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ABSTRACT 

We apply the multipole vector framework to full-sky maps derived from the first year 
WMAP data. We significantly extend our earlier work showing that the two lowest 
cosmologically interesting multipoles, t — 2 and 3, are not statistically isotropic. These 
results are compared to the findings obtained using related methods. In particular, we 
show that the planes of the quadrupolc and the octopole are unexpectedly aligned. 
Moreover, the combined quadrupole plus octopole is surprisingly aligned with the 
geometry and direction of motion of the solar system: the plane they define is per- 
pendicular to the ecliptic plane and to the plane defined by the dipole direction, and 
the ecliptic plane carefully separates stronger from weaker extrema, running within a 
couple of degrees of the null-contour between a maximum and a minimum over more 
than 120° of the sky. Even given the alignment of the quadrupole and octopole with 
each other, we find that their alignment with the ecliptic is unlikely at > 98% C.L., 
and argue that it is in fact unlikely at > 99.9% C.L. Most of the £ = 2 and 3 multipole 
vectors of the known Galactic foregrounds are located far from those of the observed 
sky, strongly suggesting that residual contamination by such foregrounds is unlikely 
to be the cause of the observed correlations. Multipole vectors, like individual aim, are 
very sensitive to sky cuts, and we demonstrate that analyses using cut skies induce 
relatively large errors, thus weakening the observed correlations but preserving their 
consistency with the full-sky results. Similarly, the analysis of COBE cut-sky maps 
shows increased errors but is consistent with WMAP full-sky results. We briefly ex- 
tend these explorations to higher multipoles, noting again anomalous deviations from 
statistical isotropy and comparing with ecliptic asymmetry manifest in the WMAP 
team's own analysis. If the correlations we observe are indeed a signal of non-cosmic 
origin, then the lack of low-^ power will very likely be exacerbated, with important 
consequences for our understanding of cosmology on large scales. 
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1 INTRODUCTION 

Following a series of increasingly successful experiments 
measuring the anisotropy of the cosmic microwave back- 
ground (CMB), full-sky maps obtained by the Wilkinson 
Microwave Anisotropy Probe (WMAP) in its first year 
of ob servation have revolut ion ized the study of the CMB 
sky jBennett et alJ l2003allhl iHinshaw et'al] [iooli)- In 
particular, a number of cosmological parameter s hav e 
been determined with high accuracy (Spergol ct al] |2003l) . 
Moreover, a WMAP-type survey opens a unique window 
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to the physics of the early universe jPeiris et alJ l2003h 
and enables tests of the standard inflationary picture, 
which predicts a CMB temperature anisotropy pattern 
that is nearly scale-free, statistically isotropic and, to 
the accuracy of all current or planned CMB experi- 
ments, Gaussian random (terms that we shall define 
more carefully below). Consequent ly, WMAP da t a led 
to m a ny s tu dies of Gaus si anity iKomatsu et al.| l2003l: 
I Parld l2004l: IChiaiig et al. l2003t iM aguciio fc Mcdciroj 
boo i IVielva et alJ l20oi iMukheriee fc Wane. i200l 
McE wen et al.l I2OO5I: ICruz et al.1 l2005l: lEriksen et alj 
[2005: 'Larson fc Wandclt '2004': 'Naselskv et al.l l2005dlbF 
iToiciro ct a l .^ ^2005: Cayon ct al. .2005) and statistical 
isotropy jEriksen et"all l2004bl: iHaiian fc Souradeed l2003l: 
iHaiian et alJl2005l :lH £Ui£u^^Sourad^eDll2005nHans^^ 
l2004allbt iPrunet et alJ 120051: iDonoehue fc Donoehud I2OO5I) 
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of the CMB. While there was no significant evidence for 
the violation of Gaussianity and statistical isotropy in the 
pre-WMAP era, some of the aforementioned post-WMAP 
studies found evidence for violation of either Gaussianity 
or statistical isotropy, or both. 

Arguably, the biggest surprises were to be found in 
temperature anisotropies on the largest angular scales. 
Most prominent among the "low-£ anomalies" is the near 
vanishing of the two-point angular correlation function 
C{d) at angular separat ions greater than about 60 de- 
grees (jSpergel et aLllioO.^f l. confirming what was first mea- 
sured using the Cosmic Background Explorer's Differen- 
tial Microwave Radiometer (COBE-DMR) a decade ago 
(Hinshaw et al. .li)96 ). Beyond this long-standing anomaly 
in the overall amplitude of the large-angle fluctuations, it 
has also been noted that the octopole of t he CMB is pla- 
nar and oriented parallel to the qua drupole ijXegmark et alJ 
l2003l:lde Oliveira-Costa et alJbOoil . Furthermore, three of 
the four planes determined by the quadrupole and octopole 
are orthogonal to the ecliptic at a level inconsistent with 
Gaussian random, statistically isotropic skies at 99.8% C.L., 
while the normals to these planes are aligned at 99.9% 
C.L. with the direct ion of the cosmological dipole and with 
the equinoxes iSchwarz et al 2004) . These peculiar correla- 
tions presumably are connecte d to north-south asym metry 
in the angular power spectrum ("Eriksei^^l and in 

the statistics of the extrema (Wandc lt et al.ll20 04'l. Finally, 
there is a suggestion that the presence of preferred direc- 
tions in the microwave background mult ipoles extends be- 
yond t he octopole to higher multipoles jLand fc Magueiia 
120053) and that there is an associated mirror symmetry 
(Land & Maaucijo 2005c). 

The correct explanation of these unexpected correla- 
tions of the low-^ features of the microwave background 
with each other and with the solar system is currently not 
known. There are four possibilities: (1) there is a system- 
atic error (an error in the data analysis or instrument char- 
acterization), (2) the source is astrophysical (i.e. an unex- 
pected foreground), (3) it is cosmological in nature (e.g. an 
anisotropic universe, such as one with nontrivial topology), 
or (4) the observed correlations are a pure statistical fiuke. 
The evidence for many of these low-^ correlations is strong 
(> 99% C.L.) and presented with a var iety o f different meth- 
ods (s ee , for examp le, Togmark et al. (2003); Schwarz et alj 
J2004I): iLand fc Magueiia (■2005a1 : ..Hansen et al,. (..2004bF : 
lEriksen et'aLr(l2004tJ) ). It is therefore unhkely that all of 
them are mere accidents. In this paper, we will at least 
attempt to clarify which of these and other correlations is 
likely to be a statistical fiuke. Whatever the source of these 
correlations, until they are understood, cosmological infer- 
ences drawn from low-^ WMAP data (including polarization 
data), such as the possibility of early reionization, should be 
viewed with a healthy dose of skepticism. 

This paper has two principal goals. First, we would like 
to further develop the theory of multipole vectors, a new 
represent ation of microwav e background anisotropy intro- 
duced bv lCopi et al.l (l2004) . In particular, we derive several 
new results, and comment on various ideas and results on 
the tests of Gaussianity and statistical isotropy that were 
subsequently presented by others using the multipole vec- 
tors or their variants. Our second major goal is to present a 



detailed analysis of t he large-angle correlations discussed in 
ISchwarz et al.l ll2004l) and extend them in several directions. 

The paper is organized as follows. In Sec. |21 we intro- 
duce the WMAP maps and their treatment used in our anal- 
ysis. In Sec. 13 we review the theory behind the multipole 
vectors, derive several new results, and comment on other 
recent work on this topic. Section|l]deals in more detail with 
the morphology of the quadrupole and octopole. In Sec. |H| 
we introduce various statistics to quantify the low-^ cor- 
relations, extend our analysis to higher I, and present the 
results. Section |S| includes a detailed analysis of the issue 
of foregrounds and Sec. |7] describes the cut-sky reconstruc- 
tion algorithm and results based on reconstructed maps. In 
Sec.|H| we present a comparison to COBE. The correlations 
of the WMAP angular power spectrum with the ecliptic 
plane is discussed in Sec. |51 We conclude in Sec. ^1 



2 THE MICROWAVE BACKGROUND SKY 

The temperature fluctuations of the microwave background 
AT{9,(j);x,t) are, in principle, functions of both the direc- 
tion of observation {0,(j)), and the location, x, and time, t, 
of the observation. In practice, essentially all our observa- 
tions occur within the solar system and over a cosmologi- 
cally short period of time, so we can ignore any local spatial 
and temporal variations in the microwave background. Our 
observations in each wavelength band are therefore the in- 
tensity of the microwave background radiation as a function 
of direction on the celestial sphere. 

2.1 The standard representation 

The most common representation of a real scalar function, 
f{0, (j)), on the sphere is as an expansion in terms of multi- 
pole moments 

oo 

f{e,<p) = Y.Me,<p) (1) 

where the i-th multipole, fi{9, <j)), is typically given by 

e 

fi{e,(P)= aemYi^{e,(P). (2) 

Here Yem{d,(t>) are the standard spherical harmonic func- 
tions. For fi{0,4)) real the complex coefficients aem need to 
satisfy the reality condition — (— l)'"af,_m and for fixed 
£ we have 21+1 independent, real degrees of freedom. 

Standard cosmological theory predicts that the CMB 
fluctuations sample a statistically isotropic, Gaussian ran- 
dom fleld of zero mean. "Gaussian random" means that the 
real and imaginary parts of the aem are each an independent 
random variable that is distributed according to a Gaussian 
distribution of zero mean. In principle, the variances of these 
Gaussian distributions (which, by the nature of a Gaussian, 
fully characterize the distribution) could be different for each 
£ and m and for both the real and imaginary parts. "Statisti- 
cally isotropic" means that, instead, these variances depend 
only on £. Thus, the expectation of any pair of aim is 

{Cl-e'm'<^lm) ~ CeSe'l^rn'm, (3) 
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where Ce is the expected power in the ^th multipole and its 
(standard) observable estimator is 

m=-l 

The set of estimators {Ce \ £ — 0, ... , 00} is caUed the angu- 
lar power spectrum. The cosmic variance in the estimators 
is 

Yai{Ce) = 2C!/{2e. + l). (5) 

Since the distributions of Gaussian variables are com- 
pletely determined by their means and variances (and since 
the aim have zero means), if the microwave background sky 
is indeed a realization of a Gaussian random, statistically 
isotropic process, then all of the accessible information in 
a microwave background temperature map about the un- 
derlying physics is contained in the angular power spec- 
trum. (Non-linear growth of fluctuations will cause depar- 
tures from Gaussianity. However, these departures are small 
at the large angular scales we will be considering.) 

As mentioned in the introduction, many studies have 
been done looking for evidence of non-Gaussianity or de- 
viations from statistical isotropy in the microwave back- 
ground. A significant difficulty is that in the absence of par- 
ticular models the range of possible manifestations of non- 
Gaussianity and statistical-anisotropy is enormous. Statis- 
tics that measure one manifestation well, can be entirely in- 
sensitive to another. Limits on particular statistics should be 
viewed in that light. Another conceptual and practical diffi- 
culty is separating tests of Gaussianity from tests of statis- 
tical isotropy. If statistical isotropy, (ap^raim) oc SiieS^i^, 
is violated, then it is unclear how one can test Gaussianity, 
since each aim could have its own independent distribution 
for which we are provided just one sample. 

2.2 Full-sky maps 

In this work, we use three full-sky maps based on the orig- 
inal single-frequency WMAP maps. The first two, the In- 
ternal Linear Combination (ILC) map and the Lagrange 
Internal Linear Combination (LILC) map, are minimum- 
variance maps obtained from WMAP's single- frequency 
ma ps bv iBennett et ai] ll2003all (the WMAP team) and 
bv lErikser^^iLr(l2004ah . r espectively. Th e third map is 
the cleaned full-sky map of iTeemark et all i2003l) (hence- 
forth the TOH map). The full-sky maps may have residual 
foreground contamination that is mainly due to imperfect 
subtraction of the Galactic signal. Further more, the full- 
sky m aps have complicated noise properties jBennett et alJ 
l2003a^ that make them less than ideal for cosmological tests. 
While one can, in principle, straightforwardly compute the 
true (full-sky) multipole vectors from the single-frequency 
maps with the sky cut (as explicitly shown for an isolatitude 
sky cut in Sec.|7J, a Galaxy cut larger than a few degrees 
will introduce a significant uncertainty in the reconstructed 
multipole vectors and consequently any statistics. 

Nevertheless, there are compelling reasons why we be- 
lieve the results of the analysis of the full-sky maps. First, all 
of the results are robust with respect to choice of the map, 
despite the fact that the ILC and LILC maps were obtained 
using a different method than the TOH map. This suggests 



that the full-sky maps indeed resemble the true Galaxy- 
subtracted microwave sky. Furthermore, one can show that 
the dominant component of the bias inherent in creating 
an ILC-type map has a quadrupolar pattern, and in par- 
ticular looks like the spherical harmonic Y2Q with an am- 
plitude of < 20 fiK (H.K. Eriksen, private communica- 
tion). For the quadrupole and octopole the dominant con- 
tributions to the Galactic foregrounds are a linear combi- 
nation of 120 and Re(y3i), which is easily seen from the 
symmetry of the Galaxy (approximately north-south and 
east-west symmetric with a hot spot at the Galactic cen- 
ter and a cold spot at the anti-center). Indeed, the sum of 
the synchrotron, free-f ree and dust WMAP foreground maps 
iBennett et alJl2003a^ in the V-band gives a2o ~ —217 fiK 
and Re(a3i) = 88 ^K. These two modes make up about 90% 
of the power in the quadrupole and octopole foregrounds. 
The next important modes in the WMAP foreground maps 
turn out to be Re(l22) and Re(V33). The quoted 20 fiK un- 
certainty in the ILC-type map for the quadrupole can thus 
be understood as a 10% uncertainty in the understanding of 
our galaxy. The uncertainty in t he octopole is conside rably 
smaller (< 10 /iK). As shown in lSchwarz et al.l ll2004f) and 
further demonstrated in Sec.|S] this Galactic contamination, 
even if present, would lead to Galactic and not the observed 
ecliptic (Solar System) correlations. 



2.3 Methodological differences between the 
minimal variance maps 

The guiding principle in the construction of the WMAP ILC, 
TOH and LILC full-sky maps is to search for a tempera- 
ture map with minimal variance. For the convenience of the 
reader we summarize the essential steps that lead to the 
cleaned full-sky maps (see the original papers for full de- 
tails) . 

In all three approaches, the input are the five WMAP 
frequency maps (K, Ka, Q, V and W band). The bands 
differ in a number of ways including noise properties and 
angular resolution. While the ILC and LILC map reduce all 
five bands to the K band resolution, the TOH map makes 
use of the higher resolution of the higher frequency bands. 
Each map can be written as 

r(i/i) = TcMB + Trcsidual(!^i)- (6) 

A combined map is created by the linear combination 

5 

T = TcMB + WiTrcsidualil^i), (7) 

with = 1. Now it is assumed that the residuals (noise 

and foregrounds) and the CMB are uncorrelated, so that 

Var(T) = Var(TcMB) + Var w;iT,esiduai(i^»)^ • (8) 

The idea is now to determine the 4 independent weights by 
minimizing the variance with respect to all pixels of a region 
of the sky. The ILC and LILC maps use the same 12 regions 
of the sky, whereas the TOH map uses 9 regions. For the ILC 
and LILC maps the weights are constants within a region 
of the sky, whereas for the TOH map the weights depend 
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on multipole number I (the minimization is done in spheri- 
cal harmonic space instead of pixel space). The regions are 
selected according to the level of foreground contamination. 
To produce the final maps a Gaussian smoothing is applied 
to soften the edges. We see that the TOH map on the one 
hand and the ILC and LILC maps on the other hand use dif- 
ferent procedures. The ILC and LILC maps differ mainly in 
the detailed implementation of the method. The drawback 
of the minimal variance method is that large CMB fluctu- 
ations tend to be canceled by large artificial foregrounds in 
order to obtain a small variance. 

We shall see that, despite the differences between the 
three maps, the results of our analysis are very similar. This 
robustness argues strongly in favor of the legitimacy of using 
these full-sky maps for multipole vector analysis. 

2.4 Kinetic quadrupole correction 

By far the largest signal in the microwave background 
anisotropy is the dip ole, recently measured by WMAP 
l|Bennett et alJl2003bjl to be (3.346 ± 0.017) mK in the di- 
rection {I = 263? 85 ± 0?1,6 = 48?25 ± 0?04) in Galactic 
coordinates. This is nearly two orders of magnitude larger 
than the root-mean-square (rms) anisotropy in the dipole- 
subtracted sky, and so thought not to be of cosmological 
origin, but rather to be caused by the motion of the so- 
lar system with re spect to the rest frame def ined by the 
CMB. As shown bv lPeebles fc WilkinsonI ljl968h . the dipole 
induced by a velocity v is T{v/c) cos6, where 9 is measured 
from the direction o f motion. Given f = (2.725 ± 0.002) K 
jMather et alJIlQQSt) . one infers that v ~ 370 kms"^ 

The solar motion also implies the presence of 

a kincmatically induced Dopple r quadrupole (DQ ; 

[Peebles & Wilkinson 1 96S ; Kamionkowski fc Kno:3 l2003l) . 
To second order in /3 = w/c ~ 10"'^, the specific intensity 
of the CMB for an observer moving with respect to the 
CMB rest-frame includes the usual monopole term with a 
black-body spectrum (oc x^/[e^ — 1], where x — hiy/lkBT]); 
a dipolar term oc cos6', linear in /3, with a dipole spectrum 
(cx x''e^/[e^ ~ the same as for primordial anisotropies) ; 
and a quadrupolar term oc 3cos^ ^ — 1, quadratic in /3, with 
a quadrupole spectrum (oc x^e'^fe^ -I- l]/[e^ — 1]^). Higher 
multipoles are induced only at higher order in /3 and so can 
be neglected. 

To first approximation the quadrupole spectrum dif- 
fers very little from the dipole spectrum across the fre- 
quency range probed by WMAP. The DQ is itself a small 
contribution to the quadrupole. It has a total band-power 
of only 3.6 /xK^ co mpared to 123.4 /xK^ from the cut- 
sky WM AP analvsis llHinshaw et alJl2003all . 195.1 /xK^ ex- 
tracted ('de Oliveira-C ost^^^m20Mr ^mn the ILC map 
(Bcrmctt et al. 2003b), 201.6 ^J.K'^ from the TOH map 
(TcEmarketalj 12003)) or 350.6 fiK^ from the LILC map 
lEriksen ^^Lril2004ai'l ) . Therefore, it is a good approxi- 
mation to treat the Doppler-quadrupole as having a dipole 
spectrum plus a small spectral distortion which we shall ig- 
nore. We can then readily subtract the DQ from any mi- 
crowave background map. 

The kinetic quadrupole is a very small contribution to 
the total theoretical power in the quadrupole, however, it 
gives rise to non-negligible contributions to some of the 
a2m- This is due partially to the low power in the observed 



quadrupole and partially to the "orthogonality" of the cor- 
rection (the correction is m dependent and the largest cor- 
rections are not to the largest a2m)- This is a well known, 
well understood physical correction to the quadrupole that is 
often ignored. This correction must be applied when study- 
ing the alignment of the quadrupole, leaving it out intro- 
duces a correctable systematic error. Though it has little 
effect on the power in the quadrupole, it has a noticeable ef- 
fect on the quadrupole orientation as shown below. We again 
stress that it is the orientation of the correction that makes 
it important. A quadrupole correction of this size pointing 
in a random direction would typically not lead to noticeable 
alignment changes. 

Figures Q |21 and El show the Doppler corrected 
quadrupole, the octopole, and their sum, for the TOH map 
in Galactic coordinates in MoUweide projection. The corre- 
sponding multipoles from the ILC and LILC maps are very 
similar. 



3 MULTIPOLE VECTORS 

Multipole vectors provide an alternative representation of 
the complete set of information represented by the full-sky 
map. They represent a radical reorganization of the infor- 
mation that nevertheless retains the integrity of the indi- 
vidual multipoles. Moreover, the information is arranged in 
such away as to be independent of the choice of coordinate 
system. Therefore the multipole vectors may be a superior 
representation for probing the null hypothesis of statisti- 
cal isotropy, and for looking for signatures of specific effects 
that might pick out special directions on the sky due to 
nonstandard inflationary physics, systematic artifacts in the 
map-making, unexpected foregrounds, cosmic topology, de- 
viations from General Relativity, or other unknown effects. 

3.1 Definition 

In the multipole vector representation, fe{9,<l)) is writ- 
ten in terms of a scalar, A^^' and £ unit vectors, 



{i, 



He A) 



n 



(9) 



Here e is the (radial) unit vector in the [6, (j)) direction. 
(Henceforth, we will use e and {d, <j)) interchangeably.) In 
Cartesian coordinates, e — (sin 9 cos <j), sin sin (j), cos 8) . Te 
is the sum of all possible traces of the first term; render- 
ing the full expression traceless. In this context, a trace 
means replacing a product of dipoles (f)'-^''-' • e) (■i)^*'-'' • e), 
by • ■y^^'^-'y Equation (|UJ can also be written as 



■ e . 

H-i ■■■1/1 



(10) 



where 



denote the trace free tensor product, and the 



sum over repeated indices is assumed. These unit vectors 
fj(id) g^j-g only defined up to a sign (and are thus "headless 
vectors"), as a change in sign of the vector can always be 
absorbed by the scalar A^^\ We have chosen the convention 
that multipole vectors point in the northern Galactic hemi- 
sphere, although when plotting them we often instead show 
the southern counterpart for clarity. 
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Note that only A^^^ depends on the total power Ce- The 
multipole vectors are independent of Ct — if all the aim of 
a given £ are multiplied by a common factor, then A too 
will be multiplied by that factor and the will remain 

unchanged. In particular, let 

atm = aem/VCe; (11) 

then, the multipole vectors, depend only on the dim 

and not on Ce. This is true independent of any assumptions 
about Gaussianity and statistical isotropy. Note, however, 
that Ce and A'*' do not contain identical information; Ce 
is the two-point correlation function while A'*' contains the 
two-point correlation as well as a particular combination of 
the higher order moments. An equivalent definition of the 
estimator Q over the full sky is 



Ce 



1 



2e + i 



dn [fe{e,(b)f 



In terms of the multipole vectors (|5J , we find 



{2e+l)Ce^[A^'^]' J 



n 



e-Te 



dQ. 



(12) 



(13) 



Evaluating this for the monopole, dipole and quadrupole we 
obtain 



Co = 47r[A(")]^ 



3Ci 



47r 



1(1)1 



3 L J ' 
5C. = %[A^^r[^ + U^ 



.(2,1) , .(2,2)^ 



(14) 
(15) 
(16) 



Notice the presence of the term ■y'^'^' • w^^'^' in the 
quadrupole expression. This implies that A'^' cannot be 
extracted solely from the angular power spectrum, but re- 
quires additional information, such as higher order correla- 
tion functions. 

A simple alg orithm for constru cting these vectors has 
been provided bv lCopi et al.l ll2004ll which builds upon the 
standard spherical harmonic decomposition. The algorithm 
relies on the observation that a dipole defines a direction 
in space, that is, a vector, and, in general, the ^-th multi- 
pole is a rank £, symmetric, traceless tensor. This tensor can 
be written as the symmetric trace-free product of a vector 
and a symmetric trace-free rank £ — 1 tensor. This proce- 
dure can be repeated recursively. It leads to sets of coupled 
quadratic equations for the components of the vectors and 
the remaining tensor whic h can be so l ved n umerically. The 
details of this are given in lCopi et alJ (|20o3). We have used 
the freely available implementation of this algorithm^ for 
our work here. 

Recently the multipole vector representation has been 
studied and used in various ways. The multipole vectors 
can be understood in th e context of harmon ic polynomi- 
als llKatz fc Weeks! 120041: iLachieze-RevI l2004f) . which has 
led to an alternative algorithm for determining the com- 
ponents of the vectors as roots of a polynomial. Expres- 
sions for A''-point correlation functions of these vectors 
(for Gaussian random aem) have been derived analytically 



1 See http:/ /www. phys. cwru.edu/projects/m pvectors/| for access 
to the code and other information. 



by iDennij (l2005h . Finally the application of these multi- 
pole vectors to the cosmic microwave background is be- 
ing actively pursued (see for examp le tCoBi_et_al 



Schwarz et al.N2004l):ISlosar fc SeliakI ll2004l):IWee'k3 Im 



2004 



Land fc Magueiid ^l2005d^ : lBielewicz et al.U2005D '). 

The multipole decomposition includes all the available 
information. Some representations related to the multipole 
vectors do not fully encode this information as discussed 
below. 



3.2 Maxwell multipole vectors 

James Clerk Maxwell, in his study of the properties of the 
sphe rical harmonics , introduced his own vector representa- 
tion (lMaxwelill89lh . In this representation, the l-th multi- 
pole of a function is written as 



MO A) = 



A 



^,(<,i) 



(17) 



where the unit vectors t)' 
tipole vectors. It is well 
unique (see lDennisl „2004 
shown that th ese vectors 
constructed bv lCopi et al 
Maxwell's representation 
techniques (see appendix 



are known as the Maxwell mul- 
known that this representation is 
and references therein). It can be 
a re pre cisely the multipole vectors 
] i2004ft . To see this we can rewrite 
usin g standard Fo urier integration 
A of lDennijIiool as 



e)+B 



(18) 



where A^^'> = (-1)'^ [(2^ - 1)!!] and B^^' is an object 
with maximum angular momentum 1—2 that ensures the 
traceless nature of fi. That is, S''^' = -A'-'^'^Ti. (Note that 
the Maxwell multipole vector representation is manifestly 
symmetric in the t)'*'"'^) This should b e compared to th e 
discussion in Sec. Ill and Appendix A of lCopi et al.U2004D . 
Thus, the constru ction outlined above and given in detail in 
ICopi et alj (|20o3) is actually an algorithm for quickly and 
easily converting a standard spherical harmonic decomposi- 
tion into a Maxwell multipole vector decomposition 11811 . 



3.3 Relation to "angular momentum dispersion" 
axes 

The peculiar alignme nt of the quadrupole and octop ole was 
first pointed out by Ide Oliveira-Costa et al.l (|20o3). They 
identified rie as the axis, which, when chosen as the fiducial 
0-axis of the coordinate system, maximizes the "angular mo- 
mentum dispersion" of each multipole. 



(AL)i 



2| 

m \ai„ 



(19) 



By varying the direction of the fiducial 2:-axis over the sky, 
and recomputing for each such choice the aem and (AL)|, 
they were able to find the choice of z-axis that maximizes 
(AL)|. This angular momentum axis distills from each mul- 
tipole a limited amount of information, reducing the 2£ de- 
grees of freedom i n the (t)^^'-'-'} to just two. 

Interestingly. Ide Oliveira-Costa et alJ (1200411 found that 
n2-n-j, = 0.9838 in the TOH map. If the £ = 2 modes and the 
£ = 3 modes are indeed statistically independent, then this 
degree of alignment has only a 1.6% chance of happening 
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Figure 1. The £ = 2 multipole from the iTegmark et alJ 1200.^) cleaned map, presented in Galactic coordinates, after correcting for the 
kinetic quadrupole. The solid line is the ecliptic plane and the dashed line is the supergalactic plane. The directions of the equinoxes 
(EQX), dipole due to our motion through the Universe, north and south ecliptic poles (NEP and SEP) and north and south supergalactic 
poles (NSGP and SSGP) are shown. The multipole vectors are plotted as the solid red (dark gray in gray scale version) symbols for 
each map, ILC (circles), TOH (diamonds), and LILC (squares). The open symbols of the same shapes are for the normal vector for each 
map. The dotted line is the great circle connecting the two multipole vectors for this map. The minimum and maximum temperature 
locations in this multipole are shown as the white stars. The direction that maximizes the angular momentum dispersion of any of the 
maps coincides with the respective normal vector as discussed in the text. 



0.019 



accidentally. (Arguably this probability should be increased 
by a factor of two, to 3.2%, since we would have been equally 
surprised if n2 had been nearly orthogonal to na.) 

We can calculate n2 explicitly for the quadrupole in 
terms of the multipole vectors. The quadrupole function is 



r^t \ l-i^.l) ~\l-{2,2) ■.\ 1~(2,1) ~{2,2) 

Q{r) = [y^ ■ ' ■ ej [v^ ' ' ■ ej ' ■ ' ' , 



(20) 



where e = r/r. We now apply to Q{r) the angular momen- 
tum operator 

L = -ire x V, (21) 
finding 

- iLQ = (-e'^'i) X e)(i)(2.2) . e) + ({)(2>2) ^ g)(^(2,i) . (22) 
We want to maximize 



(AL)l = I \n- LQ\^de 



(23) 



sky 



over all possible unit vectors to find n2. The quantity (n • 
LQ) is easily calculated: 

ifi LQ = e • [(*('•'' X n){v^^-^^ ■ e) 

+ (i)('^^' X n)(w<^''' • e)] . (24) 

The integral in 1231 is straightforward in any basis. For an 



= (sin X cos i/j, sin x sin tp, cos x) • 



(25) 



in a coordinate system where is identified with the x- 

axis and is taken to define the a;j/- plane with ' ' ■ 



(AL) 



15 L 



' X (2 + sin^ uj + cos(2[i; — 2?/;) 



+ cos(2i/'))l . 



(26) 



The partial derivative of (AL)^ with respect to x van- 
ishes at X = 0, 7r/2 and vr. We find that when x ~ 7i"/2 there 
are four minima of {AL)^ for ^ = iu/2 + mT/2, n = 0, 1, 2, 3. 
The directions x = and n are maxima. Since ip is not 
defined when sinx = 0, these are zeroes of all directional 
derivatives. Thus, the "maximum angular dispersion" is ob- 
tained in the direction normal to the plane that is de- 
fined by the two multipole vectors of the quadrupole, i.e. 



±{v 



(2,1) 



r,{2.2))/|^{2,l) 



The minima are defined 



by the directions ±{v'-^-'^^ ± v'-^-^^)/y^. 

We had previously identified the "area vectors" 



X 



(27) 

and the corresponding normalized directions to'^'*'^' as phe- 
nomenologically inte resting. I ndeed, most of the interestingj 
statist ical results of ICopi et al. ( 2004) and ISchwarz et alJ 
(2004) relate to statistics of the dot-products of area vec- 
tors (or their normalized versions) with one another or with 
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physical directions on the sky, rather than to the statis- 
tics of the dot-products of the multipole vectors themselves. 
We see that for £ = 2, the maximum angular momen- 
tum dispersion (MAMD) axis is parallel to the area vector, 
i.e. nz = ±u;(2'i'2\ 

This relation between the area vectors and the MAMD 
axis cannot extend precisely to higher £. An octopole, for 
example, has three multipole vectors |i = 1,2,3} 

which define three distinct planes with area vectors 



-,(3,1) 



1, 2, 3; i < j}. The octopole function is 



-,(3,3) 



(28) 



where the trace term is 

% = ie . . ■&(3.2)){,(3,3) + (^(3,2) . i,(3,3))^(3.1) 



(^(3,3).^(3.1)^^(3,2)^ 



(29) 



Once again, the function (n ■ LO) is easily calculated for 
some arbitrary axis n 

in-LO = e- {[(w(=''^).e)(«(=''2).e)«(^'^) 

_i({,(3.1).«(3,2))^(3,3) 



-I- cyclic permutations] x n} . 



(30) 



The MAMD axis, h^, is obtained by integrating the square 
modulus of this function over the full-sky de, and maxi- 
mizing with respect to the choice of axis n. Unlike in the 
quadrupole case, we don't find ng to be a simple combina- 
tion of the three octopole vectors v^''^''\ Rather, ria is some 
combination of the three directions which is difficult to de- 
termine analytically. 

In the specific case where the three octopole vectors lie 
in a plane, (corresponding to a so-called "planar" octopole) 
the cross product of any two are orthogonal to all three. For 
a planar octopole, the MAMD axis and the area vectors of 
the three octopole planes are all parallel. In the case of the 
observed microwave background, the octopole is relatively 
planar, for example for the TOH map; 



I (3;1,2) , (3;2,3) , (3;3,1)| 



0.8 (Iw 



(3;l,2)| 



+ \w 



(3;2,3)| 



+ \w 



(3;3,1), 



(31) 



Moreover, the three octopole area vectors surround the 
quadrupole area vector (see Fig.|HJl, so that the MAMD axis 
— which is some average of the area vectors — is even closer 
to the quadrupole axis. For these reasons, the MAMD axis is 
a moderately good representation of the octopole area vec- 
tors. This explain s that the alignment of the qua drupole and 
octo pole seen bvlde Oliveira-Costa et al] ||2004) and by us 
iSchwarz ct a^]|2004^ are indeed the same effect. 



3.4 Relation to minima/maxima directions 

The directions toward the minima and maxima of each mul- 
tipole are simple to see and at fir st glance appear to b e 
an equally useful representation (see lGluck fc Pisandl20o3) . 
This, however, is not the case. 

Consider again the quadrupole 11201 . The extrema of Q 
on the sphere occur where the angular momentum is zero. 
That is, they are the solutions of 

= -iLQ(r) 



= (t)*^''' X e){v 



(2,2) 



) + {v^^'^'^ X e){v 



(2,1) 



(32) 



The solutions of these equations occur (by inspection) along 
the directions 



-,(2,±) 



-(2,1) , --,(2,2) 



^(2,1)^^(2,2) 



Iv'-'-^^ ±V^'-^^\ ^2(l±i)(2>l).*{2>2)) 



(33) 



It is easily seen that ■i)'^'"'"' ■ t)'^'"^ = 0. Thus we see that 
for the quadrupole the maxima and minima are orthogonal 
to each other and lie in the same plane as the multipole 
vectors. The quadrupole multipole vectors bracket the two 
hot spots — the maxima occur half way between the two 
multipole vectors. 

The quadrupole multipole vectors v^-^'^^ and t)*^^'^-* con- 
tain 4 pieces of information and thus fully specify the shape 
of the quadrupole (the fifth piece of information being the 
amplitude A^'^>). The same is not true of the v^^'^\ the 
directions to the maxima and minima. Since the minima 
and maxima are orthogonal, these vectors contain only three 
pieces of information — for example, the direction of the first 
maximum, and the orientation of the plane of the maxima 
and minima. The ratio of the amplitudes of the two max- 
ima and minima is unspecified. That is, we do not know the 
relative strengths of the maxima and minima. 

For higher multipoles the minima/maxima directions, 
unlike the multipole vectors, persist in not containing the 
full information. The minima and maxima directions will 
again lack information about the relative amplitude of the 
extrema. There are also additional problems concerning the 
definition of the minima/maxima directions: In the case 
of a pure Yeo mode, there are degenerate rings of min- 
ima/maxima, which do not allow to assign a unique direc- 
tion. Moreover, the number of minima and maxima for a 
fixed value of £ is not unique. For a pure mode (spherical 
harmonic) the number of extrema depends on £ and m in 
general. To see that it is instructive to remember that £ cor- 
responds to the total number of nodal lines and m counts 
the number of meridians that are nodal lines. With this rule 
in mind one can easily see that, e.g. Y32 has four maxima 
and four minima, whereas Yjja has three maxima and three 
minima. For a general multipole there is no rule on how 
many minima and maxima we should expect and thus the 
minimum/maximum directions are only of limited use. Con- 
versely, multipole vectors always contain 2£ pieces of infor- 
mation, including information on the location, number and 
relative amplitudes of the minima and maxima. 

Therefore, while one might have at first imagined the 
minima/maxima directions to be independent, they are ac- 
tually strongly correlated. This weakens the statistical power 
of tests of the distribution of these directions. For these rea- 
sons, the statistical properties of the minima/maxima direc- 
tions are not considered further in this work. 



3.5 Relation to Land-Magueijo vectors 

iMaeueiiol il995l) discussed an alternative approach to the £- 
th multipole, which is well known to be a representation of a 
symmetric trace-free tensor of rank £. Much as for the multi- 
pole vectors, in this approach one recasts the £-ih multipole 
asa3x3x---x3(i? dimensional) Cartesian tensor, O'l"'* 
(with in = 1,2,3). One then realizes that the information 
this Cartesian tensor encodes can be recast as 2£ — 2 scalars. 
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Figure 2. The £ = 3 multipole from the iTegmark et al.l i200.'j) cleaned map, presented in Galactic coordinates. The solid line is the 
ecliptic plane and the dashed line is the supergalactic plane. The directions of the equinoxes (EQX), dipole due to our motion through 
the Universe, north and south ecliptic poles (NEP and SEP) and north and south supergalactic poles (NSGP and SSGP) are shown. 
The multipole vectors are the solid magenta (medium gray in gray scale version) symbols for each map, ILC (circles), TOH (diamonds), 
and LILC (squares). The open symbols of the same shapes are for the normal vectors for each map. The dotted lines are the great circles 
connecting each pair of multipole vectors for this map. The light gray stars are particular sums of the multipole vectors which are very 
close to the temperature minima and maxima of the multipole. The solid black star shows the direction of the vector that appears in the 
trace of the octopole, T3 J^, of the TOH map. The solid magenta (again medium gray in the gray scale version) star is the direction 
to the maximum angular momentum dispersion for the octopole, again for the TOH map. 



the 21 — 2 independent invariant contractions of the rank I 
trace-free symmetric tensor, and three degrees of freedom 
associated with an orthonormal frame. (The orientation of 
the 2-axis of the frame is two degrees of freedom. The orien- 
tation of the X-axis within the plane orthogonal to the z-axis 
is the third. The y-axis is then fixed by orthonormality and 
the convention of right-handedness.) 



iLand fc Magueiid (l2005b^ discussed the two scalars 
and the frame associated with the microwave background 
quadrupole. The two scalars are the power-spectrum and 
the bispectrum. The vectors of the frame are the eigen- 
vectors of the 3x3 Cartesian quadrupole tensor. This has 
the advantage that, at least for the quadrupole, one has 
clearly separated the issue of non-Gaussianity (the bispec- 
trum) from that of statistical anisotropy (the frame) . Land 
and Magueijo claimed that the frame they f ound was inde- 
pende nt of the multipole vectors described in'Schw arz et alJ 
i2004l) . However, Schwarz and S tarkman (private commu - 
nication, 2004 as referenced in iLand fc Magueiid bOOSdh 
showed that, in fact, two of the axes of the frame are parallel 
to the sum and difference of the quadrupole multipole vec- 
tors v^^'^' ± The third is, of course, parallel to their 
cross-product. 

Explicitly, for a quadrupole with multipole vectors 



■y'^'^-* and ■fi'^'^', the Cartesian tensor representation of the 
quadrupole 11201 is 



(34) 



.(2.1). (2,2) (2,1). (2.2) 

— 5,, 

This 3x3 matrix can be diagonalized. The three eigenvectors 
making up the Land and Magueijo orthonormal frame are 



(LM) 


-(2.1) 


X 




-(2.1) 


+ ^,(2.2) 




-(2.1) 


_i,(2.2) 



(35) 



It is noteworthy that Wj'^^' and 1)3'"'^' are in the di- 
rections of the quadrupole maxima (or minima), while 
^(LM) ^-^^ quad r upole a rea vector w^'^'^'^K (The claim by 
iLand fc Magueiid ll2005bft that these vectors did not coin- 
ci de with the quadrup ole multipole vectors or area vector 
of lSchwarz et alJ ll2004l) arose because ISchwarz et alJ (120041) 
corrected for the kinetic quadrupole as discussed in Sec. 12. 41 . 

A difficulty that arises with the Land and Magueijo 
approach is how to go beyond the quadrupole. 
iLand fc Magueiid ll2005dl) note that the 2.^ — 2 inde- 
pendent scalars associated with the ^-th multipole are just 
the 21 — 2 linearly independent combinations of dot-products 
of the I multipole vectors. The orthonormal frame is then 
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Figure 3. The i = 2 + 3 multipoles from the lTegmark et alJ l200,'jl cleaned map, presented in Galactic coordinates. This is a combination 
of Figs. and 121 with only the multipole vectors for the TOH map shown for clarity. The solid line is the ecliptic plane and the dashed 
line is the supergalactic plane. The directions of the equinoxes (EQX), dipole due to our motion through the Universe, north and south 
ecliptic poles (NEP and SEP) and north and south supergalactic poles (NSGP and SSGP) are shown. The £ = 2 multipole vectors are 
plotted as the solid red (dark gray in gray scale version) diamond and their normal is the open red (dark gray in the gray scale version) 
diamond. The £ = 3 multipole vectors are the solid magenta (medium gray in gray scale version) diamonds and their three normals are 
the open magenta (medium gray in the gray scale version) diamonds. The dotted lines are the great circles connecting the multipole 
vectors for this map (one for the quadrupole vectors and three for the octopole vectors). The minimum and maximum temperature 
locations of the 1 = 2 multipole are shown as the white stars. The light gray stars are particular sums of the £ = 3 multipole vectors 
which are very close to the temperature minima and maxima of the octopole. The solid black star shows the direction of the vector that 
appears in the trace of the octopole, T3 I29i . The solid magenta (again medium gray in the gray scale version) star is the direction to 
the maximum angular momentum dispersion for the octopole, again for the TOH map. 



constructed out of the multipole vectors, as they were for 
the quadrupole. However, there are £ multipole vectors and 
so l{£ — 1) different ways to construct the frame (not to 
mention a much larger number of ways to choose the set of 
independent scalars). To define a unique orthonormal frame 
[Land & Maguciio (2005e) chose the two vectors that have 
the two largest values of 

Ki = J2i*'-*jf- (36) 

This is not a unique choice and it is unclear how to fairly 
choose such a frame. When they claim that it is this or- 
thonormal frame that one should use to test statistical 
isotropy, we must ask which of the man y possible allowed 
choices of frame is the correct one. As iLand &i Maeueiid 
J2005el) point out the choice of frame is crucial and their 
scheme is increasingly sensitive to small fluctuations in the 
positions of the multipole vectors as £ grows. This leads 
to discontinuous noise in the Euler angles of the orthog- 
onal frame making it dif ficult to interpret the results. As 
iLand fc MagueiiolT2005(J) also point out different orderings 
of the vectors will be sensitive to different features in the 



data. Without a unique and well understood prescription 
this approach does not lead to further understanding. The 
m ultipole vectors, on the other hand are unique. It is true, 
as ICopi et all i2004) and ISchwarz et al.l lj 2004) point out, 
and as reiterated by iLand fc^aguenor(l2005dll , that they 
contain information on both statistical isotropy and non- 
Gaussianity. Optimal separation of these two properties re- 
mains an open problem. 

3.6 Polydipoles 

The directions of the multipole vectors are somewhat dif- 
ficult to interpret physically. We would argue that in large 
measure this is due to the removal of the traces done in con- 
structing the multipole vector expansion @ . An alternative 
approach would be to write the general function f{8, <p) in 
a slightly different way: 

00 00 ^ 

f{e, <p) = Y. = E 4" n^p''"' • (37) 

e=o 1=0 1=1 

The Ve's, which we term polydipoles, are similar to the 
multipoles fe in equation Q , except for our failure to remove 
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the traces in their construction. But it is this failure which 
makes them easily visualized: VeiO, 4>) is the simple product 
of £ dipoles. Therefore there are £ great circles on which 
'Pi{9,(f)) vanishes — the great circles are normal to each of 
the polydipole vectors p'^''-*. These are clearly related to the 
£ nodal curves of Yim- 

Unfortunately the polydipole expansion 1371 is unstable 
in the following important sense. Suppose that in expand- 
ing f{9,(l}) we at first neglect all contributions with angu- 
lar momentum greater than L and calculate ft and Ve for 
£ = 1, . . . ,L. If we then increase L, the fe's do not change 
(since the spherical harmonics are an orthogonal basis) but 
the Ve generally do. Thus, the polydipole expansion is well- 
defined only if the angular power spectrum falls sufficiently 
quickly at large £. The precise condition for convergence of 
the polydipole expansion has not been determined. 

The "leading" term in a multipole fe (by which we mean 
the part excluding Te in is a polydipole, which we may 
term the associated polydipole of that multipole. (Note that 
the multipole vectors are not however the vectors one would 
get in the polydipole expansion, so the £ih polydipole in a 
polydipole expansion is not the associated polydipole of the 
£th. multipole in a multipole expansion.) To the extent that 
it is easier to visualize the polydipole, it is of some limited 
interest. 



4 THE STRANGE PROPERTIES OF THE 
QUADRUPOLE AND THE OCTOPOLE 

As we have already remarked, the microwave background 
anisotropies on large angular scales seem to have several un- 
usual properties. Most widely known is that the power in the 
quadrupole, C2, is substantially less than is expected from 
the models that fit the rest of the angular power spectrum 
(and other) data. The power in the octopole is also less than 
expected, though within cosmic variance error bars. By £ = 4 
the power in the CMB is entirely consistent with theoretical 
expec tations. This was first found by COBE (.Bennett et al . 
Il996fl and has now been confirmed by WMAP. The precise 
statistic al significance of this deviation is a matter of some 
dispute llEfstathioull2004l:ISlosar et alJl2004 feielewicz etHil 
l2004l:IO'Dwver et alJbOolT 

It has also been known since COBE, b ut had largely 
been forgotten until confirmed by WMAP dSoergel et alJ 
|20o3), that the two-point angular correlation function of the 
microwave background 

c{e) = (r(ei)r(e2)) (38) 

(where ei • 62 = cos 9) is nearly zero at angu lar scales be- 
tween about 60° and 170°.* Spergel et aP i2003l) argued that, 
given the best fitting ACDM model, this is unlikely at the 
0.15% level. What has been under-appreciated, is that this 
vanishing of C{6) is not merely due to the lack of quadrupole 
power, but also due to the lack of octopole power, and m aybe 
even to the ratio of C2 : C3 : C4 iLuminet et al.ll2003^ . 

These anomalies relate exclusively to the power in the 
various multipoles. More recently, attention has turned to 
the "shapes" , "phase relationships" or "orientations" of the 
multipoles, i.e. to the information contained in the multi- 
pole vectors. As remarked above, several groups of authors 
have noticed particular anomalies. In this section, we will 



Table 1. Multipole vectors, and oriented area vectors, 

for the quadrupole and octopole in Galactic coordinates 
{I, b). All vectors are given for the TOH cleaned map after correct- 
ing for the kinetic quadrupole and the coordinates are consistent 
with how they are plotted in Figs. mOl andlSl The magnitudes for 
the oriented area vectors are also given. 



Vector 


I 


b 


Magnitude 


^,(2,1) 


118?9 


25? 1 




i,{2.2) 


11?2 


16? 6 




i„(2;l,2) 


74? 3 


-56? 6 


0.990 


i,(3,l) 


86? 9 


39? 3 




i,{3,2) 


22? 6 


9? 2 




£,(3,3) 


-44?9 


8? 2 






101?6 


-49? 8 


0.902 


it,(3;2,3) 


-6?3 


-79? 5 


0.918 


it,(3;3,i) 


38? 4 


-38? 9 


0.907 



first describe them, mostly in the language of multipole vec- 
tors, before proceeding to try to assign them some statistical 
significance in the next section. 

We believe that the observed lack of power at large an- 
gles justifies singling out the two most responsible multipoles 
for particular scrutiny. To do so is no more dubious than the 
practice of firefighters to respond to the house where the fire 
alarm is ringing rather than to all the houses in the neigh- 
borhood. That these lowest multipoles represent a physically 
interesting scale — the recent scale of the horizon, especially 
the scale of the horizon at approximately dark energy dom- 
ination — makes it doubly justified to focus our attention, 
at least initially, on them. 

Table contains the multipole vectors and area vec- 
tors for the quadrupole and octopole that will be discussed 
below. 

4.1 The queerness of the quadrupole 

Figure shows the £ = 2 multipole from the lTeemark et all 
(2003) cleaned map after subtraction of the kinetic 
quadrupole. The solid line is the ecliptic plane and the 
dashed line is the supergalactic plane. The directions of the 
equinoxes (EQX), dipole due to our motion through the 
Universe, north and south ecliptic poles (NEP and SEP) 
and north and south supergalactic poles (NSGP and SSGP) 
are shown. The multipole vectors are plotted as the solid 
red symbols for each map (see figure caption), while the 
open symbols of the same shapes are for the normal vector 
for each map. The dotted line is the great circle connecting 
the two multipole vectors for this map. The minimum and 
maximum locations of the temperature in this multipole are 
shown as the white stars. The following observations can be 
made about the quadrupole: 

(i) The great circle defined by the quadrupole multipole 
vectors ?)'^'^' and f)'^'^' passes through the NEP and SEP. 
This is especially true for the TOH and LILC maps, with 
some slight deviation for the ILC map. We can rephrase this 
to say that the (normalized) quadrupole area vector u;^^'^'^-' 
lies on the ecliptic plane. 

(ii) The axis of this great circle (i.e., w^'^'^'^^) is aligned 
with both the dipole and the equinoxes. 
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4.2 The oddness of the octopole 

Figure 13 shows the ^ = 3 muhipole from the lTegmark et all 
l|200 ,T) cleaned map. Many of the features are the same as 
for the quadrupole map (see Fig. Ql. Note that there are 
now 3 multipole vectors (closed symbols) and 3 normal vec- 
tors (open symbols) plotted. Also there are 3 dotted lines 
showing the great circles connecting each pair of multipole 
vectors for this map. Several perplexing observations can 
also be made about the properties of the octopole. 

(i) The octopole has three multipole vectors: 

^,(3.1)^ ^(3.2) 

and v'^^'^K One of these, v'^^'^\ lies quite near the ecliptic 
poles. Therefore, the two great circles defined by the pairs 
(^)(3,i)^ ^(3,2)^ and the pair (u^^'^', -S^^'^') each nearly pass 
through the ecliptic poles. But in fact, one of these great 
circles passes much closer to the poles than the position 
of t)''''^'. The associated area vectors, ly''^'^'^' and ty'^'^'^', 
therefore lie on (or nearly on) the ecliptic. 

(ii) The third pair of multipole vectors, [v'^^''^\v^^'^^) de- 
fine a great circle that includes the supergalactic poles. This 
pair also lies within 9° of the Galactic plane. The associated 
area vector tu'^'^'"^' , therefore lies on the supergalactic plane 
and only 10?5 from the Galactic poles. 

(iii) The octopole is noticeably planar, but not over- 
whelmingly so. Namely, the area vectors of the octopole 
cluster only somewhat compared to expectations. 

(iv) The three maxima and three minima are very nearly 
at ±v^^'^^ ±t)'^'^^ ±v^^'^\ even though these are not generi- 
cally where the angular momentum is zero. (Note that there 
are 8 such choices but only 6 are either minima or maxima 
direction.) This contrasts with the quadrupole, where the 
extrema can be shown analytically to be at 

(see section . 

(v) The vector from the octopole trace, Ta, (see Ean. l29ll 
lies on the ecliptic plane. The importance of this vector is 
unclear. Note that Ta is proportional to the difference be- 
tween the octopole and the polydipole, "Pa. 



4.3 The remarkable relation of the quadrupole 
and the octopole 

Figure shows the sum of £ = 2 and £ = 3 in the TOH 
cleaned map, after kinetic quadrupole subtraction. We can 
see that, beyond the separate oddities of £ = 2 and £ = 3, 
there are additional unexpected relationships between the 
quadrupole and octopole: 

(i) The quadrupole normal vector i2;'^'^'^^ (open red sym- 
bols) — or equivalently the quadrupole MAMD axis n2 — 
is aligned with the octopole MAMD axis ns (magenta star). 
It is also aligned with the three octopole area vectors (open 
magenta symbols), especially with one of them. 

(ii) The ecliptic carefully traces a zero of the combined 
map, doing so almost perfectly over the entire hemisphere 
centered on the Galactic center, and even relatively well over 
the antipodal hemisphere. 

(iii) Two of the extrema south of the ecliptic are clearly 
stronger than any north of the ecliptic. The weakest south- 
ern extremum is essentially equal in power to the strongest 
northern extremum. (It is slightly stronger in the LILC map, 
slightly weaker in the TOH and ILC maps.) 



Table 2. Values of various vector dot-products for the TOH 
cleaned map, corrected for the kinetic quadrupole. We show the 
values of the dot-products of the three octopole area vectors with 
the quadrupole area vector (A^, i = 1,2,3); the dot-products of 
the three octopole normal vectors with the quadrupole normal 
vector (Di, i = 1,2,3; the normals are the unit area vectors); 
and the dot-products of the four quadrupole and octopole area 
and normal vectors with the north ecliptic pole (NEP), the north 
Galactic pole (NGP), the north supergalactic polo (NSGP), the 
dipole and the equinox. 



Value of the dot-product 



Tost 


i„(2;i,2) 


i„(3;i,2) 


iu(3;l,3) 


iu(3;2,3) 


A, 






0.851 


0.783 


0.762 


w ■ 


NEP 


0.027 


0.161 


0.041 


0.481 


w ■ 


NGP 


0.827 


0.688 


0.570 


0.903 


w ■ 


NSGP 


0.392 


0.262 


0.630 


0.0011 


w ■ 


dipole 


0.974 


0.883 


0.755 


0.674 


w ■ 


equinox 


0.968 


0.886 


0.681 


0.766 


D, 






0.953 


0.872 


0.838 


iv ■ 


NEP 


0.027 


0.179 


0.045 


0.523 


w ■ 


NGP 


0.835 


0.763 


0.629 


0.983 


iv ■ 


NSGP 


0.396 


0.291 


0.694 


0.0012 


iv ■ 


dipole 


0.984 


0.979 


0.832 


0.733 


iv ■ 


equinox 


0.978 


0.982 


0.751 


0.834 



5 STATISTICS 

We now revisit the statistics used to quantify the align- 
ment of the vari ous vectors w ith each other and with the 
Solar System. In ISchwarz et ah. I|2004), we considered the 
dot-products, Ai, of the three octopole area vectors with 
the quadrupole area vector, and the mutual dot-products, 
Di, of the three octopole normal vectors (area vectors nor- 
malized to unit length) with the quadrupole normal vector. 
We found that both the Ai and the Di are unusually large. 
Note that we can take as a convention that Ai ^ since 
each vector is ambiguous in sign. The variou s dot-products 
are sh own in Table |21 As we emphasized in ISchwarz et alJ 
(^2004), all of the aforementioned alignments are statistically 
significant at 99.8% C.L. or higher. 

To compute all probabilities we compare statistics ap- 
plied to WMAP maps to that applied to Monte Carlo sim- 
ulations. Unless otherwise noted, Monte Carlo simulations 
are comprised of 100,000 realizations of Gaussian random, 
statistically isotropic maps with WMAP's (inhomogeneous) 
pixel noise. 

5.1 Requirements for robust statistics 

In previous work (|Schwarz et al.l2004l . we used a particular 
set of statistics, for example, for a set of vectors that have 
dot-products {Ai \ i = 1, . . . ,n} (with one another, or with 
a particular physical direction on the sky) with "unusually" 
high values. We asked what is the probability that a ran- 
dom Monte Carlo map has the highest dot-product higher 
than Ai, the second highest one higher than A2, the third 
highest one higher than A3 , and so on down to the n-th such 
dot product. For several cases that number was, as reported, 
very small. For example, for the dot products between the 
quadrupole area vector w^^'^'^^ and each of the three oc- 
topole area vectors only 21 out 
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of 100,000 MC maps (for the TOH DQ-corrected map) sat- 
isfied the criterion, i.e. only 21 had Ai larger than the TOH 
value of Ai, A2 larger than the TOH value of A2 and A3 
larger than the TOH value of ^3. 

One can ask if this statistic preferentially returns a 
small probability even if there is really none to be found. 
This is clearly not the case — the vector dot-products could 
have been either large, small or "average"; we noticed they 
were large and found an easy to understand statistic that 
quantified the effect. Weeks (private communication) how- 
ever has pointed out that the above-described statistics for 
the Ai do not define an ordering relation on the set of pos- 
sible Ai; they therefore implicitly incorporate some a pos- 
teriori knowledge. One would therefore like to confirm this 
result with different, independent statistics. 

5.2 S and T statistics — definitions 

To quantify the various alignments we found, it is desirable 
to choose the statistics in such a way that the a posteriori 
knowledge of the particular nature of the alignments is not 
used to find unjustly small probabilities. With that in mind 
we define and discuss two statistics, S and T, which do define 
ordering rel ations. The first of th ese was briefiy mentioned 
and used in lSchwarz et all (1200411 as suggested by Weeks. 

Two natural choices of statistics which define ordering 
relations on the three dot-products Ai, each lying in the 
interval [0, 1], are: 

S = ^{Ai+Ai + As) (39) 

and 

T = 1- 1 [(l-Ai)V (1-^2)2 + (l-yla)'] . (40) 

Both 5* and T can be viewed as the suitably defined "dis- 
tance" to the vertex {Ai,A2, A3) — (0, 0, 0). A third obvious 
choice, {Ai + Ai + Al)/S, is just 25" - T. Of course many 
other choices exist, involving higher powers of Ai. 

One could also ask about the probability that, for ex- 
ample, two out of three normals are aligned, and so we gen- 
eralize the definitions to 

m 

g{n,m) _ ly^^^ (41) 
i=l 

m 

m 

1=1 

where the Ai are ordered from largest to smallest Note that, 
if the Ai are large (near 1), both S'"'™' and T'"'™' will be 
large (near 1). 

One could further generalize these statistics to arbitrary 
weighting by different alignments; for example 

n / ^ 

^(n.o) _ ^^^^ / (43) 

1 = 1 / 2 = 1 

with ^ ^ 1. Hereafter we consider only the values 
Qi = or 1. Finally, there is nothing special about the area 
vector products Ai and we apply the statistics S''"'™-' to 
dot-products of the vector normals, Di, and also to dot- 
products of the normals with specific directions or planes in 
the sky (ecliptic plane, NGP, supergalactic plane, dipole and 



equinox) that were discussed in lSchwarz et al.l ll2004ll . When 
comparing to planes, we order the dot products (taken with 
the plane's axis), from smallest to largest. 

5.3 S and T statistics — quadrupole and octopole 

Figure |1] shows histograms of S*"'""* statistics (left column) 
and statistics (right column) for Gaussian random, 

statistically isotropic MC maps, as well as values for the 
TOH DQ-corrected map. The statistics shown are 

5(3,™) 

and t'^'™' for the intrinsic alignment of the octopole area 
vectors with the quadrupole area vector, and S*-'*'™-' and 
alignment of normals with the ecliptic plane. 
Galactic poles, and supergalactic plane. The dipole and 
equinox alignments are similar to the ecliptic plane and are 
not shown. 

Note several interesting features in Fig. 0| First, the 
probabilities for the algebraically-related values of 5^"'^^ 
and r("'i' are by definition identical, since they measure 
the extremeness of a single parameter. Second, evidence for 
alignment of the quadrupole and octopole and for alignment 
of multipoles with the ecliptic plane (and similarly with the 
dipole and equinoxes) is strongest when all dot-products are 
considered — i.e. when m = n. This gives us further confi- 
dence that these probabilities are generally not dominated 
by one or two unusual alignments, but rather when align- 
ment of all four normals, either mutual or with the specified 
direction, are considered. 

We note that the alignment of normals with either the 
supergalactic plane or the Galactic poles is dominated by 
a single normal, w'^^''^'^^ . In the case of the supergalactic 
plane, this normal is only 0?07 away from the plane (in the 
TOH map), while the other three are not particularl y close 
to the plane at all; see Table |5| Figure and Schwa rz et alJ 
i2004). Therefore S^^'^' < 5''''''") for m ^ 2. Although this 
normal is still within approximately 1?5 in the ILC and LILC 
maps, this is still sufficient to raise the probability of S'*'^' 
to approximately 10%. The fact that only the S''''^' statistic 
and only the TOH map show small probabilities suggests 
that the supergalactic correlation is a statistical fiuctuation. 

The same normal uj''''^'^' is approximately 10° from the 
Galactic pole. This is largely responsible for the measured 
correlation of the quadrupole and octopole with the Galactic 
pole, although it is true that 5^'*''*' < 5'*'^' for this case. We 
discuss this further in Sec. 15.41 

In the absence of any model, 5"'"'"' and T^"'") statis- 
tics seem the fairest choice, and the one we adopt. These 
statistics treat planes and directions identically (i.e. the or- 
dering of the relevant dot-products does not matter). Table 
l^shows the probabilities for the seven different tests applied 
on the TOH, LILC and ILC maps using the statistic. 
We show the results based on both DQ-corrected and DQ- 
uncorrected maps. 

Finally, there has been some confusion regarding con- 
straints from the application of the S statistic to the nor- 
malized, Di, and unnormalized, Ai, area vectors. It has been 
claimed that the S statistic applied to the area vectors Ai is 
both not a measure of alignment between the octopole and 
quadrupole and not a very robust statistic. It was argued 
that one should only use the Di. The heuristic argument is 
that S applied to the Ai is not actually a measure of align- 
ment since the Ai incorporate information about the lengths 
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Figure 4. Histograms of 5'"'™' statistics (left column) and T^"'™) statistics (right column) for Gaussian random, statistically isotropic 
Monte Carlo maps with m = 1 . . . n of the most aligned area vectors considered separately in each panel. First row panels show the 
n = 3 mutual dot-products of quadrupole and octopole area vectors Ai , second row shows dot-products of the n = 4 normals with the 
ecliptic plane, third row shows dot-products of the n = 4 normals with the north Galactic pole while the fourth row shows dot-products 
of the n = 4 normals with the supergalactic plane. Specific values for the WMAP (TOH DQ-corrected map), for each m, are shown with 
vertical lines. The numbers show the percentage of MC maps that have a more extreme value (i.e. larger mutual dot- products Ai, larger 
products with the NGP, or smaller dot-products with the NEP and NSGP). In other words, the numbers show the extremeness of each 
vertical line's value in the corresponding histogram. Note that the statistical significance is strongest when all vectors are considered 
(that is, when m = n), except for the supergalactic plane where a single octopole normal is only 0?07 away from this plane while the 
other three are not particularly unusual. 
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Table 3. Table of probabilities (in percent) for the seven different tests using the TOH, LILC and fLC maps. We show the results based 
on both DQ-corrected and uncorrected maps. The statistics 5'"'"', computed for each microwave background map and compared to 
100,000 realizations of Gaussian random, statistically isotropic maps with WMAP's pixel noise, has been used to compute all probabilities. 
Except for the test Di, all tests are based on the dot-products of the area vectors, w, of the quadrupole and octopole. 



Test 


TOH DQ-corr 


LILC DQ-corr 


ILC DQ-corr 


TOH uncorr 


LILC uncorr 


ILC uncorr 


A, 


0.117 


0.602 


0.289 


0.582 


2.622 


0.713 




1.246 


1.309 


2.240 


1.262 


1.309 


2.567 


ecliptic plane 


1.425 


1.480 


2.006 


1.228 


1.735 


2.724 


NGP 


0.734 


0.940 


0.508 


0.909 


1.265 


0.497 


SO plane 


14.4 


13.4 


8.9 


11.6 


10.2 


6.5 


dipole 


0.045 


0.214 


0.110 


0.093 


0.431 


0.207 


equinox 


0.031 


0.167 


0.055 


0.064 


0.315 


0.080 



of t he area vecto rs, not just their directions (see, for exam- 
ple, |^eki|20oi). However, what this means is that the S 
statistic applied to the Ai weights the contribution of each 
plane according to how well the two associated multipole 
vectors define that plane — the closer they are to orthogo- 
nal, the more well defined the plane is, and the more heavily 
the plane is weighted; the more nearly parallel they are, the 
less well defined the plane is and the smaller its weighting 
in the statistic. This weighting thus seems entirely intuitive 
and appropriate. 

The confusion noted above has arisen from inconsis- 
tently applying the test to maps, some with the kinetic 
quadrupole correction applied and others without it. The 
specific concern about the robustness of this statistical test 
was that the S statistic applied to the Ai of the LILC map 
was apparently much less unlikely than the S s tatistic ap- 
plied to the Ai of the TOH map JWeeksll20oil . However, 
this result compared the Doppler-corrected TOH map to the 
uncorrected LILC map. The former has probability around 
0.1%, whereas the latter has probability 2.6%. When we 
Doppler-correct the LILC map (as one should), we find that 
the probability falls to around 0.6%, comparable to that of 
the TOH map (see Tablej^l. Also the Doppler corrected ILC 
map yields a similar 0.2% probability. Far from challenging 
the robustness of the S statistic for the area vectors, the 
probabilities derived from different maps support it. 

5.4 The evidence for ecliptic (and other) 
alignment 

The values in Table |21 show alignment of the area vectors 
of the quadrupole and octopole with each other and with 
some specific physical directions. Also shown in Table El is 
the statistical evidence for a correlation with the dipole and 
equinox direction (at larger 99.7%C.L. in all three maps), 
with the NGP (at larger 99%C.L.) and the ecliptic plane 
(at larger 98%C.L.). As discussed in Sec. 15.31 our use of the 
5*'''*'''' statistic (instead of the S^^'^-* statistic) shows that the 
correlation with the supergalactic plane is not significant. 

Previously we have claimed that there is evidence for 
ecliptic alignment, but not for Galactic alignment. Yet Ta- 
ble 121 shows a slightly higher significance for the Galactic 
than the ecliptic alignment. In this Section we demonstrate 
why the ecliptic correlation is significant and the Galactic 
one is not. Simultane ously, we consider the suggestion by 
iBielewicz et al.l (|200^ that the only important correlation 
is between the quadrupole and octopole area vectors them- 



selves. We show that actually there is a > 98% C.L. fur- 
ther alignment with the ecliptic plane (and argue that this 
alignment is in fact > 99.9% C.L. unlikely), but that the 
additional alignment with the Galaxy is not significant. 

To begin we take the correlations in the quadrupole 
and octopole (their area vectors) to be fixed as measured. 
We then compute 100,000 rotations of the TOH quadrupole 
and octopole on the sky and compute the S'^'^''^^ statistic 
for each direction in each of these rotate d microwave back- 
ground skies (see IBielewicz et al.l ll2005l) for a similar ap- 
proach). The results are shown in Fig. |K1 The histogram is 
the distribution of the S'-^'*-' we get from these TOH rotated 
skies and the dashed vertical lines are the values of S^^'^'^ for 
each of the ecliptic plane, NGP, supergalactic plane, dipole 
and equinox axes. In Table |31 we list the percentiles of the 
values of s'^'^''^^ for these five physical directions for all three 
of the TOH, the ILC and the LILC maps. (The distributions 
for the ILC and LILC maps are quite similar to the TOH 
distribution.) 

The results are striking. The percentile for the ecliptic 
plane is between 0.2% (LILC) and 1.7% (ILC). The per- 
centile for the Galactic pole is between 87% and 90%, so that 
the two-sided probability is only between 74% and 80%. (We 
justify below why a two-sided probability is not appropri- 
ate for the ecliptic plane alignment.) This shows that, given 
the observed shapes and alignment of the quadrupole and 
octopole, the evidence from the area vectors for additional 
correlation with the ecliptic is at least 10 times stronger 
than for additional correlation with the Galaxy. There is 
also mild evidence for additional correlation with the dipole 
or equinox at approximately the 95% C.L. (approximately 
the 90% C.L. when we take two-sided probabilities). 

Qualitatively we can understand from inspection of 
Fig.|21why the quadrupole and octopole normals are so much 
better correlated with the ecliptic than with the Galactic 
pole. These four normals essentially surround the ecliptic, 
therefore it is relatively hard to be more correlated; similar 
is true, though to a lesser extent for the dipole and equinox 
directions. On the other hand just one of the normals comes 
somewhat close to the Galactic pole, while the other three 
are further away. The situation for the dipole and equinox 
is intermediate to these two. 

It is also very important to notice that the quadrupole 
area vector does not contain all the information about the 
quadrupole. That means the information that a zero of the 
sum of quadrupole and octopole traces the ecliptic for about 
1/3 of the sky is not contained in our s'*''*'' statistics. More- 
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over, the 5''' statistics does not make use of all the infor- 
mation contained in the octopole area (normal) vectors. To 
see this consider a somewiiat idealized sky which is dom- 
inated by in some frame. (The is nearly true for the 
£ — 2 + 3 sky.) In this case, all the 1 = 3 multipole vec- 
tors lie approximately in a single plane, and the three nor- 
mal vectors are closely aligned, and define a small circle. 
Any rotation around the axis defined by the centre of that 
small circle will give rise to the same results in our statis- 
tical tests, because this rotation merely leads to a rotation 
of the Y33 minima and maxima within the octopole plane. 
The three minima and three maxima are separated by three 
nodal lines, which are great circles. For one of those great 
circles to be the ecliptic, the axis of rotation (the direction 
defined by the three normal vectors) would have to lie on 
the ecliptic. This is what we have nearly found to be the 
case in the data. 



5000 



4000- 



3000- 



2000 



1000- 




S*^"^' statistics 



However, placing the normal vectors on the ecliptic 
plane does not itself guarantee that the ecliptic plane is one 
of the null contours between extrema. The freedom to rotate 
all the multipole vectors in their plane remains. With 60° 
between extrema, the chance of the ecliptic plane being a 
null-contour within the observed tolerance of about 3? 5 is 
about 6%. Since this rotational freedom is entirely indepen- 
dent of the alignment of the area vectors, we can multiply 
the probability of the area-vector alignment by the proba- 
bility of the rotational alignment to obtain at most 0.1% 
(ILC), and as little as 0.01% (LILC). (It is not appropriate 
to use two-sided probabilities, because if the normal vec- 
tors had been aligned with the ecliptic poles rather than the 
ecliptic plane, then the ecliptic plane could not have been 
a null contour between extrema.) The ecliptic plane only 
traces the null-contour over one third of the sky due to the 
fact that the sky isn't a pure Y33 mode. However, the eclip- 
tic plane also doesn't pass between just any two extrema 
but instead splits the weaker extrema in the north from 
the stronger extrema in the north. Thus these two effects 
approximately cancel each other leaving us with the proba- 
bilities estimated above. Though this is an estimate and a 
more detailed statistical analysis is warranted, the probabil- 
ity of the quadrupole and octopole being this aligned with 
the solar system is unlikely at greater than 99.9% C.L. 

We note that it has been suggested that one should 
reduce the significance of this discovery by some large 
number of possible "physical great circles" with which we 
could have noticed a correlation — that our focus on the 
ecliptic is purely (and by implication fa tally ) a posteriori 
(Biolcwicz ct al. 2005; Slosar & Soljak 2TO^. To the con- 
trary, we would argue that given the experiment there are 
precisely two great circles and their axes with which one 
must look for correlations — the ecliptic and the Galactic 
equator. The former because WMAP orbits the sun deep 
within the solar system, and correlations with the eclip- 
tic could be a sign of either systematic errors or a local 
foreground; the latter because the Galaxy is an important 
foreground source, and correlations to the Galactic equa- 
tor would be a sign of residual Galaxy contamination (cf. 
Sec.inj. 



Figure 5. Histogram of the S^*'^^ statistics applied to the TOH 
map quadrupole and octopole area vectors and a fixed direction 
or plane on the sky, where the area vectors have been rotated 
together in a random direction 10^ times. Vertical lines show the 
S statistics of the actual area vectors applied to the ecliptic plane, 
NGP, supergalactic plane, dipole and equinox directions (Table 
El shows the actual product percentile ranks among the random 
rotations for all three full-sky maps). This Figure and Table Bl 
show that, even given the relative location of the quadrupole- 
octopole area vectors (i.e. their mutual alignment), the ecliptic 
plane, dipole and equinox alignments are unlikely at the > 95% 
C.L. while the NGP and supergalactic plane alignments are not. 



Table 4. Percentile ranks of the quadrupole-octopole area vector 
dot-products with specific directions, among the 10^ such prod- 
ucts where the area vectors have been rotated together in a ran- 
dom direction. We show dot-products for the TOH, LILC and 
ILC area vectors, and the NEP, NGP, NSGP, dipole and equinox 
directions. This Table shows that, even given the relative loca- 
tion of the quadrupole-octopole area vectors (i.e. their mutual 
alignment), the ecliptic plane, dipole and equinox alignments are 
unlikely at the > 95% C.L. while the NGP and supergalactic ones 
are not. 



Test TOH DQ-corr LILC DQ-corr ILC DQ-corr 



ecliptic plane 1.0 0.2 1.7 

NGP 87 88 90 

SG plane 34 33 25 

dipole 95.6 93.8 94.5 

equinox 96.1 94.4 96.4 



5.5 Single multipole alignment test for ^ ^ 50 

So far we have devoted attention to alignments of the area 
vectors defined by the quadrupole and octopole. We would 
now like to investigate alignment in higher multipoles. To 
simplify the analysis, we consider a single multipole £ at a. 
time. We then ask whether the multipole in question is un- 
usually planar. Note however that this test is by no means 
exhaustive in finding unusual correlations — for example, 
we are not considering pairs of multipoles, as we did for 
the quadrupole and octopole. In fact, the octopole, when 
considered alone, is not unusually planar as first shown by 
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iTegmark et al.l i2003h and further illustrated using our tests 
below. Nevertheless, considering the more complicated tests 
applied to higher multipoles is not efficient at this time with- 
out a specific suspected theory or model in mind. Moreover, 
we axe curious to investigate whether the "bites" observed 
on the angular power spectrum at a handful of multipoles 
jHinshaw et a l. 2003a) are somehow correlated to the mul- 
tipole vector alignments in those multipoles. Therefore, we 
proceed with testing the single multipole alignment of area 
vectors. 

For a fixed multipole £ we have A'^ = (.{£ — l)/2 area 
(headless) vectors Wi. We would like to find the plane that 
is defined by these vectors. We define the plane as the one 
whose normal, n, has the largest dot product with the sum 
of the area vectors. Since Wi ■ n is defined only up to a 
sign, Wi is headless, we take the absolute value of each dot 
product. Therefore, we need to find n that maximizes 



1 " 



Ni 



(44) 



An exact procedure to find the best aligned direction is non- 
trivial due to the presence of the absolute value. Instead of 
writing an absolute value we replace Wi by SiWi, where Si is 
±1, such that there are only positive terms in the sum. The 
problem is that we cannot know the signs before we know 
the best aligned directions and thus we have to try out all 
the possibilities. There are A'^^ area vectors, and thus there 
are 2^* possible choices for the set {si}. As is easily seen, for 
large £ this becomes computationally more expensive than a 
numerical search for the best aligned direction. For any fixed 
set of signs {si}, the best aligned direction n expressed in 
spherical coordinates, {6,<j}), is given by 



tan (p = 
tan 6 = 



(45) 



(TOi a; COS (p + Wi^y Sm < 



where Wi^a; denotes the x component of Wi, etc. 

In practice, we used an iterative solution to the prob- 
lem: pick some trial direction n so that we can immediately 
compute S. Repeat many times for n randomly chosen on 
a unit sphere, and simply find one that has maximal S. We 
find that the search converges after a few thousand trials for 
the direction n. 

We compute 5wmap for any given multipole £ in the 
TOH map, and then compare it to the value of at 
that multipole, computed in a large number of Monte Carlo 
realizations of a Gaussian random, statistically isotropic 
sky with WMAP's pixel noise added. We then rank-order 
iSwMAP among the S'^^q. A low rank indicates that the 
WMAP areas are roughly aligned, or equivalently, that the 
multipole vectors are planar and most of the power lies in 
the corresponding plane. Conversely, a high rank would im- 
ply that the multipole in question is non-planar relative to 
a Gaussian random, statistically isotropic expectation. We 
use 10,000 Monte Carlo realizations at Z 20 and at any / 
that has a very low or high rank in order to obtain sufficient 
statistics; at all other £ we use 1000 realizations. 

TableOshows the ranks of the S statistic for 3 < £ ^ 50, 



Table 5. Ranks of the S statistics, R, given as percentiles. Low 
ranks indicate multipoles that are planar (i.e. their area vectors 
are aligned), while high ranks indicate multipoles that are non- 
planar. The statistic 5, defined in Ean. l44l is applied to the TOH 
cleaned map and compared with Gaussian random, statistically 
isotropic skies to obtain the ranks. We used 10,000 realizations 
of such skies for £ ^ 20 and for £ = 44, and 1000 realizations for 
all other multipoles. Values within 5 percentiles of or 100 are 
shown in bold. 

Single-£ Ranks of the S Statistic (%) 
£ R£R£R£R£R 



2+3 


0.35 


11 


42 


21 


32 


31 


33 


41 


26 


2 




12 


19 


22 


48 


32 


46 


42 


3 


3 


7 


13 


3 


23 


20 


33 


65 


43 


86 


4 


75 


14 


5 


24 


15 


34 


11 


44 


99.7 


5 


99.7 


15 


62 


25 


48 


35 


13 


45 


35 


6 


4 


16 


4 


26 


46 


36 


96 


46 


31 


7 


45 


17 


3 


27 


96 


37 


89 


47 


49 


8 


77 


18 


18 


28 


80 


38 


47 


48 


27 


9 


95 


19 


19 


29 


58 


39 


86 


49 


91 


10 


65 


20 


29 


30 


4 


40 


12 


50 


19 



given as percentages^ . For example, 25% would indicate that 
25% of Monte Carlo maps are more planar and 75% are less 
planar at that £. Note that, while the octopole alone is planar 
only at the 92% C.L., the quadrupole and octopole together 
are planar at 99.65% C.L. The planarity of the quadrupole 
and octopole is therefore very significant, in agreement with 
other tests (see e.g. Ai and Di in Table ISJ. The £ = 5, in 
contrast, is non-planar and only 0.3% of Monte Carlo real- 
izations exhibited lower planarity. What has been called the 
"spheri city" of £ = 5 was first pointed out bv lEriksen et alJ 
(Hooii); their test found it unusual at the 5-10% level. 

Visual inspection of tablej^suggests that there is an ex- 
cess of both high and low values. Simple attempts to quan- 
tify this do indeed find such anomalies at between 95% C.L. 
and 99% C.L. We might therefore conclude that the align- 
ment test as defined in Eqn. I|44^ gives strong hints of some- 
thing unusual at 4 5C £ ^ 50 in the TOH cleaned map, but 
without further evidence, the case is not sufficiently strong 
to stand on its own for any bold claims. Moreover, we find 
that at £ ^ 8, the values of the 5 statistic differ substantially 
among the TOH, ILC and LILC maps. 

Apart from looking at the 5 statistics, we also inspected 
the best aligned directions n (the direction that maximize 
5). We would expect that only the best aligned directions 
of planar multipoles have a well defined meaning. For £ < 8 
only £ — 6 is singled out by our statistics. We find the corre- 
sponding vector at {I, b) = (152?4, 50?3), which is 46?2 from 
the ecliptic pole. Among the higher best aligned directions, 
£ — 21 and £ = 44 are within 2?6 and 9?0, respectively, of 
the dipole. All other vectors are more than 10° from any 
physical direction studied in this work. 

In addition to the S statistic described a bove, we also 
applied Bingham's statis tic test of isotropy (iFisher et alJ 
1993; Mo rgan et al ] l200,1 ). Let us assume we have A'^ unit 
vectors with components {xi,yi, Zi) (i = 1, . . . , N) and that 

^ The quadrupole alone has only one area vector which picks out 
a unique direction n and therefore cannot be used with the S 



statistics as defined in Eqn. 1441 . 
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we want to check whether they are distributed isotropically. 
We construct the orientation matrix 



(46) 



which is real and symmetric with unit trace, so that that 
the sum of its eigenvalues {k — 1,2,3) is unity. For an 
isotropic distribution all three eigenvalues should be equal 
to 1/3 to within s tatistical fl uctuations. Bingham's modified 
statistic B* (Bingha^^T^) is defined as 



B* 



B 



3 



157V 
2 



^^ + ^B- 
84 147 



5292 



B' 



where 



(47) 



For isotropically distributed vectors and N ^ 1, B* is dis- 
tributed as xi- Here we compare WMAP's value of B* for 
a given multipole to Monte Carlo simulations directly, and 
therefore do not require assuming A'^ ^ 1. 

Bingham's statistic results for 2 ^ ^ ^ 50 are broadly 
consistent with those for the 5 statistic shown in Table |S] 
(and hence we do not show them separately). We find that 
£ = 5 is non-planar at the 99.8% C.L. but, apart from that, 
other multipoles are neither planar nor non-planar at a level 
not expected from such a statistical sample. Finally, note 
that Bingham's statistic, being general and coordinate in- 
dependent, does not use all information and is typically not 
as strong as the coordinate dependent tests. 



5.6 Higher multipole angular momentum vectors 

The angular momentum dispersion 11911 can be maximized 
for all multipoles and thus serve as a statistic. That is we 
can find the axis ni around which (AL)^ is maximized. The 
spherical harmonics provide an irreducible representation of 
the rotation group in three dimensions and transform as 
Y'l — "FJ^D'*' where Ye is a vector of the £-th multipole 
spherical harmonics {2£ + 1 components) and D'^^ is a ro- 
tation of this multipole. Since a scalar function (such as 
the temperature) is invariant under rotations the aim must 
transform as a'^ = D^ae under rotations. The rotations can 
be parametrized in terms of the Euler angles a, /3, 7 in the 
zyz representation as 



D(a/37) 



exp(|:L.)exp(|L,)exp(^L.) 



(48) 



where Ly and Lz are the y and z components of the angular 
mo mentu m operator, respectively. The discussion here iol- 
lows lEdmonda ( ^19600. F or an alternative rep resentation see 
appendix A of Ide Oliveira-Costa et al] ||20o3)- 

To perform the maximization it is convenient to use the 
matrix representation for the rotations 



m' m 



(49) 



where 



m' m / J 



(-1)*-™'- + '^'W - -m-'W + ""OK^ - ■m)\ 



k\{l - m' - k)\{l - m - fc)!(m + m' + k)\ 



(^cos-j (sm-j 



2k-\-7n -\~m 



2£-2fc-m'-m 



(50) 



In this representation, the angular momentum disper- 
sion 11911 becomes 

(AL)?= ^ aL'awe"'™'-'""'^5]m^d(^],J/3)d^l„(/?) 

m' ,71%" m 

= y H^'] „(7)G<^ ,(/3) 

/ J m'm" ^ ' ' m" m' ^' ' 



= Tr(HW(7)G(^)(/?)) 



(51) 



Notice that this expression separates into a term that de- 
pends only on 7 and the a^m, H'*-'(7), and a term that only 
depends on /3, G^^-'(/3). To extremize this function we take 
derivatives of (AL)^ with respect to (i and 7 which also 
separates. It is easy to show that 



d.,H^^] „ (7) = i(m' - m)H^^] „ (7) 



(52) 



and that both G^',^, (/3) and dg G^f},,^! (P) can b e calcu- 
lated quickly and efficiently (see lEdmonda (Il960h for de- 
tails) . These rotation angles are related to standard Galactic 
coordinates via 



{l,b) = (7- 180°, 90° -/3). 



(53) 



Finally, we will find it convenient to work with the nor- 
malized angular momentum dispersion (the t statistic of 
Ide Oliveira-Costa et al.1 (Hooi)) 



(AL)| 



{AL)j 



(54) 



The normalized dispersion takes a value between {£+ 1)/S£ 
and one. 

We have maximized the normalized angular momentum 
dispersion for £ = 2 to 20 for the TOH cleaned map. The ILC 
and LILC maps give similar results. The results are shown 
in table |S| For each multipole we provide the direction in 
Galactic coordinates, {I, b), for the axis around with the an- 
gular momentum dispersion is maximized and the value of 
the maximum angular momentum dispersion. We have also 
performed 10,000 Monte Carlo simulations of Gaussian ran- 
dom, statistically isotropic skies and found the maximum 
angular momentum dispersion for each one of them. The fi- 
nal column in the table gives the percentage of MC skies that 
had an angular momentum dispersion larger than that for 
the WMAP data. We have also done the same procedure for 
a joint fit to the quadrupole and octopole {£ = 2-1-3). That 
is, we find the single axis that maximizes both multipoles. 

We find that the octopole is planar, but only at 89% 
C.L. The strong correlation between the quadrupole and oc- 
topole is seen by the fact that less than 0.4% of all Gaussian 
random, statistically isotropic skies have the quadrupole and 
octopole this well aligned. W e again confirm the "sp hericity" 
oi £ = 5 first pointed out bv lEriksen et al ] J2004ali . We find 
the angular momentum dispersion to be very low — all but 
0.56% of the MC skies have a value larger than the WMAP 
data. They have also suggested that £ = 6 is planar; by this 
test it is somewhat planar with only 3% of MC skies being 
more planar than the data. We find a total of 5 multipoles 
that are somewhat planar (less than 5% of MC skies having 
a larger angular momentum dispersion), those being £ — 6, 
13, 14, 16, and 17. We only find the £ = 5 muhipole to be 
particularly non-planar. 
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Table 6. The maximum angular momentum dispersion results for 
the TOH cleaned map. Shown in the table is the multipole num- 
ber, £, the Galactic coordinates of the direction in which the axis 
around which the angular momentum dispersion is maximized, 
(Z,fe), the value of the normalized angular momentum dispersion 
around this axis, (AL)'j, and the percentage of Monte Carlos that 
had a maximum angular momentum dispersion larger than the 
value found from the TOH map. Values within 5 percentiles of 
and 100 are shown in bold. See the text for details. 



"shape" statistic which finds the preferred axis and the pre- 
ferred m for each multipole. The statistic is defined as 



I 



MC Larger 



2+3 


-112?7 


59? 7 


0.962 


0.37% 


2 


105? 7 


56? 6 


0.993 




3 


121?6 


62? 


0.942 


11.24% 


4 


-106?0 


36? 3 


0.637 


75.93% 


5 


-170?0 


24? 


0.484 


99.44% 


6 


-160? 5 


44? 1 


0.817 


3.00% 


7 


-119?7 


54? 8 


0.585 


58.13% 


8 


149? 7 


20? 1 


0.546 


62.28% 


9 


162? 1 


77? 9 


0.504 


80.69% 


10 


153? 8 


11?1 


0.511 


70.17% 


11 


-76? 9 


12? 2 


0.540 


45.15% 


12 


46? 3 


37? 2 


0.565 


28.75% 


13 


-71?7 


41? 7 


0.651 


2.73% 


14 


59? 3 


9? 3 


0.636 


2.60% 


15 


-81?3 


31? 1 


0.503 


49.27% 


16 


137? 9 


78? 7 


0.632 


1.82% 


17 


-163? 5 


30? 2 


0.593 


4.70% 


18 


152? 8 


19?0 


0.540 


14.85% 


19 


-121?2 


60? 3 


0.530 


16.20% 


20 


124? 2 


20? 9 


0.500 


28.78% 



Setting aside the £ = 2 + 3 result, we see that 6 of the 
18 angular momentum dispersions are in either the top or 
bottom 5 percentile. The probability of having 6 or more of 
the 18 so anomalously high or low is 0.6%. We also see that 
of these 6, all but £ — 5 are in the top 5 percentile. The 
probability of having 5 or more of the 18 so anomalously 
high is 0.15%. 

Inspecting the directions of maximum angular momen- 
tum dispersion we find that only the £ = 4 direction is close 
to one of the physical directions under consideration: its dis- 
tance to the ecliptic pole is 10? 3. Note that this confirms the 
qualitative impressio n from looking at the £ = 4 map (see 
ISchwarz et al.1 l)2004l) ') that this mode has its minima and 
maxima aligned with the ecliptic plane. It is also interesting 
to note that the directions given in table [S] especially for 
the planar multipoles, are consistent with the ones found as 
the best aligned directions. 

These results are another suggestion that the higher £ 
multipoles are not statistically isotropic. Reassuringly, com- 
parison of Tables|S]and|H]shows the same 5 multipoles which 
had a high angular momentum dispersion also exhibited 
comparably low ranks of the S statistic, while £ = 5 showed 
a high rank of S. This difference in range of £ in Sees. 15. til 
and 15. 51 was purely a result of computational limitations. 



5.7 "Shape" statistic 

The angular momentum dispersion searches for planarity 
through a weighted averag e that favors modes with m = £. 
iLand fc Maeueiid J2005af) have suggested the use of the 



ri = maxr, 
where 



Om.,0 ) \a,£n 




(55) 



(56) 



and n is the z-axis of the coordinate system in which the a£m 
are computed. Note that the angular momentum dispersion 
is a weighted sum of these terms, 



(AL)| 



1 2 m 

m=0 



(57) 



The maximization of the shape statistic II55II follows 
the same formalism as for the angular momentum disper- 
sion and will not be discussed further (see Sec. 15.61 . We 
ha ve performed this maximi zation and confirm the results 
of iLand fc Maeueiid (|2005eJ) . In particular, we find that 
the surface defined by is complicated with many local 
maxima. The results are quite sensitive to the data and 
are not consistent among the three fu ll-sky maps made 
from t he WMAP data (see Fig. 2 of iLand fc Maeueiid 
(Hooli)). Unfortunately this sensitivity is not understood 
in terms of features of the data. That is, the variabil- 
ity in the results cannot be understood in terms of fea- 
tures such as non-Gaussianity or a violation of statistical 
isotropy. The sensitivity of the shape statistic is related to 
the difficulty in uniquely defining the Land-Magueijo vectors 
(Land & Maeuciio 20056, also see Sec. 13.51 . For these rea- 
sons, the shape statistic does not serve as a robust statistic 
for separating nor understanding Gaussianity versus statis- 
tical isotropy. 



6 FOREGROUNDS 

So far, we have taken into account the effects of noise in the 
full-sky cleaned maps by including the WMAP pixel noise 
into our Monte Carlos. We now explore the effect of the 
foregrounds on the quadrupole-octopole anomaly. 

While it has repeatedly been emphasized that there 
might be residual foreground contamination left in the 
cleaned maps, it seems that such a contamination should 
lead to Galactic and not ecliptic correlations. Here we ex- 
plicitly show this with a quantitative analysis. We slowly add 
the measured WMAP foreground contaminations to WMAP 
full-sky maps and monitor how the directions defined by the 
quadrupole and octopole change. 

Let rcMB(iT-) be the cleaned-map microwave back- 
ground temperature in some direction, and TpoR(n) the 
temperature from one of the three basic foreground maps 
(thermal dust, free-free emission or synchrotron emission) 
provided by the WMAP team. We form the total contami- 
nated map as 



Ttot(n) = TcMB(n) + cTFOR(n) 



Var(rcMB) 



(58) 



Var(rFOR) ' 

where c is the foreground fraction. Note that the second 
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Figure 6. Trajectories of the multipole vectors and tlieir normals as increasing amounts of the foreground is added to or subtracted 
from the microwave background map. We have use the synchrotron V-band foreground and the TOH cleaned map; the coefficient c 
shows the rms contribution of the foreground to the microwave background (see Eqns. and |F5J. The solid diamond symbols show 
the zero-foreground locations of the multipole vectors while the solid stars refer to their pure-foreground locations. Similarly, the open 
diamond and star symbols refer to the beginning and end points of the normal vectors. On each trajectory we label a few values of 
the coefficient c. The top panel shows 1 = 2 (two vector and one normal) and the bottom panel shows £ = 3 (three vectors and three 
normals). Not all trajectories end on plotted symbols. For these trajectories they end on the multipole vector (or normal) that is the 
negative of the plotted vector at the start of the trajectory. Note that, in the large-foreground limit, the quadrupole vectors move near 
the z-axis and the normal into the Galactic plane, while for the octopole all three normals become close to the Galactic disk at 90° from 
the Galactic center. Therefore, as expected Galactic foregrounds lead to Galactic, and not ecliptic, correlations of the quadrupole and 
octopole. 
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term has been normalized so that 

Var(Ttot) = Var(rcMB)+c' Var(rFOR)^^4^^ 

Var(JFOH,j 



= Var(rcMB)(l + c ) 



(59) 



so that the rms contribution to the total rms temperature 
from the added foreground is a factor, c, relative to the con- 
tribution of the cleaned microwave background map. (For 
reference, the constant •y^Var(TcMB)/ Var(rFOR) is of or- 
der unity in all cases we consider.) In .Copi et a l. (2004), we 
have performed tests with c ^ 0.2 and found no significant 
changes to the results in that paper. Note that the con- 
tribution in power of known foregrounds to the microwave 
background, after removing the foregrounds, is estimated to 
be less than a percent in the V and W bands (,Bennctt e t al„ 
However this estimate is for the multipole range 
2 ^ £ ^ 100 while the foreground contami nation is most 
signifi cant at low multipoles (see Fig. 10 in iBennett et alJ 
|2003^. Therefore it is reasonable to assume that the contri- 
bution of residual foregrounds at large angular scales is less 
than about 10% in power, or c < 0.3. 

We have added increasing amounts of foreground, cor- 
responding to c taking values from zero (no foreground) 
to ± 1000 (essentially pure foreground) and recomputed 
the multipole vectors and their normals. Fig. |n] shows the 
trajectories of the multipole vectors and their normals as 
increasing amount of the foreground (V-band synchrotron 
map produced by WMAP) is added to the microwave back- 
ground map. The solid diamond symbols show the zero- 
foreground locations of the multipole vectors while the solid 
stars refer to their pure- foreground locations. Similarly, the 
empty diamond and star symbols refer to the zero- and pure- 
foreground normal vectors. On each trajectory we label a 
few positive and negative values of the coefficient c. The top 
panel of Fig. |S| shows 1 = 2 (two vector and one normal) 
and the bottom panel shows £ — 3 (three vectors and three 
normals) . 

In the pure- foreground limit, the quadrupole vectors 
move near the 2:-axis and their normal into the Galactic 
plane as is expected for an almost pure Y20 mode (see the 
discussion in Sec. 12. 21 . In the same limit, two of the octopole 
multipole vectors move close to the Galactic poles and one 
close to the Galactic center. Consequently, all three normals 
become close to the Galactic disk at 90° from the Galactic 
center. This is indeed the signature of the expected ReCKsi) 
domination in the foreground octopole. Therefore, both the 
vectors and their normals clearly migrate from locations cor- 
related with the ecli ptic and other directions discussed in 
ISchwarz et al.l (|20o3) and in this paper to locations speci- 
fied by the Galactic foreground emission. Further, note that 
appreciable admixture of the foreground (|c| > 0.1 for the 
quadrupole and |c| > 0.3 for the octopole) is necessary for 
this migration to become apparent by eye. This confirms the 
argument in Sec. I2.2l that the quadrupole foreground is the 
most critical one. 

The reader might notice that two of the octopole multi- 
pole vectors are within 10° of the Galactic plane, or, equiv- 
alently, one of the octopole normals is about 10° from the 
Galactic poles. One might ask if that could be a sign of resid- 
ual Galactic contamination of the cleaned full-sky maps. 
However, as seen in figure and table 01 given the observed 



pattern of quadrupole and octopole area vectors, the align- 
ment of these area vectors with the NGP is significant at 
< 90% C.L. (compared to a 99% C.L. correlation with the 
ecliptic plane). Moreveor, since a correlation of the normal 
with the Galactic plane would have been even at least as 
noteworthy, this 90% C.L. figure should be reduced to 80%. 

We also observe that when adding foreground, one of 
the two multipole vectors close to the Galactic plane moves 
far away from that plane. The other multipole vector close 
to the Galactic centre does not move very far (about 10°). 
But, a 10° alignment of one of the cleaned-map multipole 
vectors with one of the foreground multipole vectors is not 
statistically significant by any statistical test applied in this 
work. Even less so when we consider that we would remarked 
as well on a similar alignment of a cleaned-map area vector 
with a foreground area vector. 

Finally, we do not expect more than one foreground 
multipole vector to lie near the Galactic plane since the 
dominant fairground mode of the octopole is indeed Re(y3i). 
The second biggest foreground mode is Re(l33), which has 
its three multipole vectors in the Galactic plane at I = 30°, 
90° and 150°. But, the cleaned full-sky map does not resem- 
ble that pattern either. We are not able to identify a Galac- 
tic contamination of the cleaned full-sky maps and thus see 
no evidence to question more significant ecliptic correlations 
found in the previous sections. 

We have checked that the results are qualitatively un- 
changed if we use the W-band synchrotron map, or the 
V-band free-free and dust foreground maps instead of the 
V-band synchrotron map. We have done some further test- 
ing, recomputing the S statistics applied to different align- 
ments but now with the foreground incrementally added 
to both the microwave background and Monte Carlo maps. 
We found that the results are consistent with those inferred 
from Fig.|S|and indicate that large admixture of the known 
G alactic foreground wo uld not cause the alignments found 
in lSchwarz et al.l J2004) . 

Results of this investigation are therefore in agreement 
with the intuitive expectation: Galactic foregrounds lead to 
Galactic, and neither ecliptic nor dipole, correlations of the 
quadrupole and octopole. In fact, it is difficult to see how any 
known foreground that has most of its power in the Galactic 
plane can lead to the solar system correlations that we find. 



7 EFFECTS OF CUT SKIES 



As mentioned in lCopi et al.l ll2004) . Galaxy cuts of a few de- 
grees or larger introduce significant uncertainties to the re- 
constructed full-sky multipole vectors. This is precisely why 
we used the cleaned full-sky maps — sky cuts of 20 degrees 
or so would simply lead to large uncertainties in our statis- 
tical tests. Nevertheless, we would like to look at the issue 
of sky cuts in more detail; in particular, we would like to 
explore how correlation significance varies as we introduce 
a sky cut. 

We are interested in finding the full-sky (true) aim 
which are derived from the full-sky temperature distribu- 
tion and denoted by al^- In general, we cannot view the 
full sky as the Galaxy obscures our field of view and must 
be cut out. There are well know techniques for relating the 
cut sky decomposition, al^, to the true sky decomposition 
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Figure 7. Quadrupole-octopole probabilities for the TOH DQ-corrected map for an increasingly larger isolatitude cut of ±(degrees 
shown), performed symmetrically around the Galactic plane (left column), ecliptic plane (central column) or an arbitrarily chosen plane 
(right column). We consider the S statistic probabilities applied to the ecliptic plane, north Galactic pole, dipole and the equinoxes (first 
to fourth row respectively). The solid line is the mean value, while the dark and light regions represent 68% C.L. and 95% C.L. regions, 
respectively, from 1000 realizations of reconstructed coefficients that take into account the noise in the reconstruction process. The 
dashed line denotes the probability obtained from the full-sky map, corresponding to the case of zero cut. 



JWandelt et alJl200lt iMortlock et alj|2002l) . We briefly dis- 
cuss the key facts here. The decompositions are related by 



where 



(60) 



(61) 



and Scut is the cut sphere. There are fas t, stable recursion re- 
lations for calculatin K these Wu',mm' llWandelt et aljliooil : 
iMortlock et alJl2002^ ■ In this work, we restrict ourselves to 
longitudinal cuts symmetric across the Galactic (xy) plane. 
In this case, m — m' and W is a symmetric matrix. For no- 
tational convenience we will drop the m index and keep in 
mind that the subsequent equations hold independently for 
each m. Thus we write 
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a1 = '^Wu'4,. (62) 

I' 

Since information is lost in the cut, W is not an invert- 
ible matrix. We can, however, replace W with an invertible 
matrix VV constructed from W by removing the rows and 
columns with small eigenvalues. That is, W = VAV^ where 
A is the diagonal matrix of eigenvalues such that \j = Aj 
if \j > Athreshoid and Xj = A~^ =0 otherwise. A threshold 
of Athreshoid < 0.1 is typically sufficient and is used in this 
analysis. An estimate, for the true decomposition is 

a\ — W^^}aii. (63) 
I' 

We can prevent leakage of power from non-cosmological 
monopole and dipole modes by projecting out these modes 
using a partial Householder transformation (see Appendix 
C of Mortlock ct al. 2002 for details). 

There is an error in this approximation which is evident 
from the fact that 

a* = W"^a'= =W"^Wa' (64) 

and W^^W 7^ I due to the loss of information in the cut. 
Our error in the approximation is found to be 

^|a*— d*|^^ — ^(a' — a')* (a' — a')^^ (65) 

= (I- W"'W)((a')*(a*)'^) (l-W"'W)^. 

Here {(a^)*a^) — Ci and this error can be readily calculated. 

Fig. |7| shows the quadrupole-octopole correlation prob- 
abilities with the sky cut between zero and ±20 degrees per- 
formed along the Galactic plane (left column) , ecliptic plane 
(central column) or an arbitrarily chosen plane'^ (right col- 
umn). We consider the S statistic probabilities applied to the 
ecliptic plane, north Galactic pole, dipole and the equinoxes 
(first to fourth row respectively). The solid line is the mean 
value, while the dark and light regions represent 68% C.L. 
and 95% C.L. regions, respectively, from 1000 realizations 
of reconstructed atm. coefficients that take into account the 
noise in the reconstruction process. While the increasing cut 
clearly increases error in the vector reconstruction and there- 
fore uncertainty in the final probability, it is clear that in 
essentially all cases the probabilities remain consistent with 
the full-sky values at 95% C.L., and in most cases at the 
68% C.L. for cuts up to 10 degrees. 

This figure clearly shows that sky cuts of a few degrees 
or larger introduce significant uncertainty in the extracted 
multipole vectors and their normals, leading to increased 
error in all alignment tests. Nevertheless, the cut-sky align- 
ments are consistent with their full-sky values even for rela- 
tively large cuts. Note that the shift of the mean value of the 
alignments (black curves in the panels of Fig. |7|l to less sig- 
nificant values, as the cut is increased, is entirely expected: 
an unlikely event, in the presence of noise in the data, be- 
comes less unlikely because any perturbation will shift the 
multipole and area vectors away from their aligned locations. 

^ The "arbitrary" plane is chosen as one obtained by rotating the 
map in Galactic coordinates by +45 degrees around the z-axis and 
then by —60 degrees around the new x-axis. The resulting map 
has neither the ecliptic nor the Galactic plane located along the 
equator. 



While the resul ts of this exercise are in good agreement with 
those found by ISlosar fc SeliakI (|2p04') and iBielewicz et alJ 
(2005), unlike these authors, we emphasize that the cut sky 
is always expected to lead to shift in the alignment values 
and to increased errors (see again Fig. |7|l. 



8 COMPARISONS WITH COBE 

Since the alignments we are studying are on very large 
scales (i.e. quadrupole and octopole scales), it is natural 
to ask whether the y can be seen in COBE-DMR maps 
iBennett et alJll99'^ . COBE angular resolution is about 7°, 
which is more than sufficient for this test. However, full-sky 
maps produced by the COBE team are very noisy (G. Hin- 
shaw, private communication) while, as discussed in Sec.|71 
using the cut sky maps produces too much uncertainty in 
the vectors for Galaxy cuts larger than a few degrees. 

Fortunately, Markov chain Monte Carlo (MCMC) 
Gibbs sampler realizations of the full-sky COBE maps have 
been produced . Usin g their Global Estimation Method, 
IWandelt et alj (|20^) generate 10,000 realizations of the 
COBE sky consistent with DMR measurements and ex- 
pected foregrounds. Following their own conservative ap- 
proach, we drop the first 1000 maps which might be affected 
by the burn- in of the MCMC, and choose every 200th map 
from the remaining 9000. The resulting 45 maps are essen- 
tially uncorrelated and their analysis is computationally un- 
demanding. We then compute the S statistic for each COBE 
map and rank-order it relative to 100,000 Monte Carlo real- 
izations of the Gaussian random, statistically isotropic sky. 
The COBE derived values for the statistic S is represented 
by histograms, while a WMAP full-sky map is represented 
by a single value. We then ask whether the two are consis- 
tent. Note that we have corrected all COBE maps for the 
DQ using the procedure described in Sec. 12.41 

The results are shown in Fig. |H| The a;-axis in each 
panel shows the logarithm of the probability of S (the prob- 
ability is simply its rank relative to 100,000 Monte Carlo 
realizations of the Gaussian random, statistically isotropic 
sky). The vertical lines show values for the TOH, ILC and 
LILC maps. The histogram shows values for the 45 COBE 
maps produced using the MCMC Gibbs sampler. Obviously, 
there is significant variation in the COBE statistics which 
traces to the fact that the MCMC maps are based on in- 
complete sky information. Nevertheless, we see that the sta- 
tistically significant alignments found with WMAP are in 
most cases consistent to the results from COBE. For ex- 
ample, 11 out of 45 COBE maps show the alignment with 
the north ecliptic pole that is equally or less likely than 
that in the WMAP maps, while 4 to 8 COBE maps (de- 
pending on which WMAP map is considered) show equal or 
less likely alignment with the equinoxes. Note too that the 
alignment with the supergalactic plane differs significantly 
in the three WMAP maps. Consequently, the comparison of 
MCMC-derived COBE maps with WMAP cleaned full-sky 
maps shows that COBE data are consistent with WMA P in 
regards to all alignments found in \Schwarz et al\ \200A) . 
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Figure 8. The oriented area statistics from fuII-sky WMAP maps compared to those from Markov chain Monte Carlo realizations of the 
COBE full sky based on COBE-DMR cut-sky data. We consider the S statistic applied to the dot-products of quadrupole-octopole area 
vectors {Ai) and the normals [Di), as well as the sum (i.e. the 5-statistic) of the four dot-products of the area vectors, one quadrupole 
and three octopole, with the special directions or planes — the ecliptic plane, NGP, dipole and equinox. The WMAP maps considered for 
the T OH, ILC and LILC as shown by the vertical lines. The histograms are from a total of 45 MCMC COBE maps from ; Wandclt ct al] 
J200'tfl (their maps 1000, 1200, . . . , 9800) which are sufficiently separated in the chain so as to be essentially uncorrelated. 



9 ANGULAR POWER SPECTRUM — 
ECLIPTIC PLANE VERSUS ECLIPTIC 
POLES 

There are at least three points in the binned WMAP angu- 
lar TT power spectrum that deviate from the predictions of 
the best fit ACDM model at comparable or greater statis- 
tical significance to the power deviation of the quadrupole. 
These are clearly seen in Fig. 12 of Bcmictt ot al. (2003h!) 
in the bins I = 20-24, I = 37-44 and I = 201-220, which 
are visibly low, high and low respectively. These deviations 
are approximately 2, 2.5 and 1.5 times the estimated error 
in the average Ci in each of those bins. In each case, this is 
largely cosmic variance dominated, with only a small frac- 
tion of the error being due to statistical error. Nevertheless, 
it should be appreciated that the statistical significance of 
these deviations may decline in the second or third year 
WMAP maps. 

Various explanations have been offered for these devi- 
ations in the angular power spectrum as ar i sing from fun- 
damen tal physics (e.g. iGordon fc Hul (|20o3); lEn^^^^j] 
i2005fl and references therein). However, it is instructive to 
examine the angular power spectrum computed separately 
using data from the ecliptic plane ve rsus data from the eclip - 
tic pole. This is presented in Fig. 7 of lHinshaw et al] ll2003bl) . 
[Note that this Figure has been replaced in the final, pub- 
lished version of the paper (.Hins h aw e t al. 2003a) with a 
cross-band power spectrum that shows some of the same 
features.] The figure caption instructs the reader to "note 
that some of the 'bite' features that appear in the combined 
spectrum are not robust to data excision." Specifically, the 
three deviations at £ ~ 22, 40 and 210 are found only in the 



ecliptic polar data and not in the ecliptic planar data. This 
suggests that there are some continuous parameters — the 
latitude of the planar-polar division, the orientation of the 
reference plane relative to the ecliptic, etc. — as a function 
of which the separation is approximately maximized by the 
WMAP team's choice of an ecliptic planar-polar division. 

A more detailed study is merited after future data re- 
leases. In the meantime, the dip at the first peak, while 
the least significant, has the virtue of being in a region 
of the angular power spectrum that has pr eviously been 
probed by multiple ex periments — TOCO iTorb et et all 
ll999l:lMiller et alll999l). Saskatoon jNetterfield et al.ll996l) ' 
Boor nerang (iNetterfield et alj20ollJones et al.l2005^. Max- 
ima ('Hanan^^^L*2'oOO^ . and Archeops llBenoi^'^1 JI2003I : 
iTri s tram et al. 2005). Of these, only Archeops saw evidence 
of a dip in the power at the first peak, though it was of lim- 
ited statistical significance. Other experiments, saw no such 
dip with greater statistical confidence. However, other than 
WMAP only Archeops covered a significant portion of the 
north ecliptic hemisphere; the others covered much less sky 
or the southern-sky. Thus if the dip in the first peak is lo- 
calized to the northern sky, especially to the region near the 
north ecliptic pole, all the experiments could be consistent. 
Notably, Archeops also shows evidence of excess of power at 

e ~ 40. 

Figure 7 of lHinshaw et al.l ^l2003b^ contains further per- 
plexing anomalies in the \ow-i sky — the angular power 
spectra extracted from the ecliptic planar and polar regions 
disagree for £ < 10. The differences are much larger than 
would be expected from statistics alone, these low-^ Ce's 
being very well measured (as can be seen from the error 
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bars in Fig. 8 of the same paper): 

1) There is a nearly complete absence of ecliptic polar power 
in the ^ = 6 — 7 bin that is highly reminiscent of the findings 
in tables 13 and |S| that ^ = 6 is very "planar." This suggests 
further that said planarity is closely aligned with the eclip- 
tic, and that this connection extends to ^ = 7 as well. This 
is also reflected by the fact that the best aligned vectors for 
£ = 6 and £ — 7 from Sec. 15.51 are 46° and 23°, respectively, 
from the ecliptic pole, which already indicates that there 
is little power in the (ecliptic) polar cap region defined as 
/3 > 30° in the WMAP analysis. 

2) The aplanarity of ^ = 5 (and to a lesser extent 1 = 4), 
as seen again in tables |S| and |S| may also be reflected in the 
notably even distribution of power between the plane and 
poles. 

3) There is a dramatic deficit in ecliptic planar power com- 
pared to ecliptic polar power at £ = 2 and £ = 3. 

If indeed there is contamination of the microwave back- 
ground from some source in the ecliptic (north?) polar region 
that is responsible for the deviations at £ m 22, 40 and 210, 
then it is possible that this contamination is also the source 
of the law-£ microwave background in the polar regions. In 
this case, the maximum cosmic contribution to the CMB at 
low £'s would be the lesser of the ecliptic planar and polar 
values (barring an exceptional cancellation). In particular, 
the cosmic C2 and C3 would be bounded by the ecliptic 
(planar) values. These are dramatically lower than the cus- 
tomary values as extracted from the full sky with a Galactic 
cut, and considerably less consistent with theoretical expec- 
tations than even the current low values. 



10 CONCLUSION 

The multipole vector forma l ism fi rst introduced to the study 
of the CMB bv lCopi et al.l (|20^ has proven to be a useful 
means of studying the structure of the CMB particularly 
on large scales. In this work, we have provided an extensive 
discussion of the multipole vector formalism highlighting the 
fact that the multipole vectors provide an alternative, com- 
plete representation of a scalar function on a sphere (see 
Sec. 13 for details ). In particu l ar, w e have pointed out that 
the algorithm of ICopi et alJ (|20o3) converts the standard 
spherical harmonic decomposition into the multipole vec- 
tor representation first discussed by Maxwell 11811 . We have 
shown how the multipole vector formalism relates to the 
previously studied maximum angular momentum dispersion 
(MAMD) directions, Land-Magueijo vectors, and tempera- 
ture minima/maxima directions. Note that, unlike the mul- 
tipole vectors, neither the MAMD nor the minima/maxima 
directions contain the full information of a multipole, and 
are thus not complete representations of the microwave sky. 
The Land-Magueijo vectors and scalars are a complete rep- 
resentation, but suffer from a rapid proliferation of arbitrary 
choices for £ > 2. 

As noted, the multipole vectors are an excellent way to 
study alignments and correlations in the microwave sky. We 
have provided a qualitative description of the striking prop- 
erties of the quadrupole and octopole in Sec. 2] We note that 
there are strange properties for the quadrupole and octopole 
individually as well as jointly. Not all of these unexpected 
properties are independent of each other and an explanation, 



whether statistical fluke, residual foreground contamination, 
or real CMB features, remains to be determined. 

By eye the properties of the quadrupole and octopole 
multipole vectors seen in Figs. are striking. To quantify 
these correlation we have used the S and T statistics for the 
oriented area and normal vectors (see Sec. 15.21 and applied 
them to the quadrupole and octopole (see Sec. 15. 31 . We con- 
firm the alignment of the quadrupole and octopole planes at 
greater than 99% C.L. We also confirm that the quadrupole- 
octopole planes are aligned with the geometry and direction 
of motion of the solar system. In particular, they are per- 
pendicular to the ecliptic plane at approximately 98% C.L. 
and to the dipole and equinox at > 99.8% C.L. They are 
also perpendicular to the Galactic poles at > 99% C.L. 

We have shown that the alignment with the ecliptic 
plane remains at 99% C.L. when the quadrupole-octopole 
alignment is taken as given (for the TOH-DQ map; similar 
or stronger results hold for the other maps, see Sec. 15.41 . 
The correlations with the dipole and equinox remain at ap- 
proximately the 95% C.L. However, the correlations with 
other directions, such as the Galactic poles, do not persist. 
This strongly supports the reality of the ecliptic correlation 
in the data and suggests that the aforementioned alignment 
with the Galactic poles is accidental. 

We further stress that the 99% C.L. correlation of the 
quadrupole-octopole planes with the ecliptic plane is a lower 
bound. As is evident from Fig. |3] and discussed in more de- 
tail in Sec. l5.4l the ecliptic plane carefully threads its way be- 
tween the temperature minima and maxima of the £ = 2 + 3 
map separating the weak power in the northern ecliptic sky 
from the strong power in the southern ecliptic sky. This ex- 
tra feature that is manifest in the multipole vectors is not 
contained in our statistics of oriented area (nor normal) vec- 
tors. Thus the analyses discussed in this work and in the 
literature which rely solely on dot-products of the oriented 
area (or normal) vectors are not using all the information 
available in the multipole vectors. Dot-products of oriented 
area (and normal) vectors are well suited for identifying and 
defining planes but do so at the expense of the informa- 
tion of the structure with respect to these planes. (For the 
quadrupole the area vector contains only three out of four 
pieces of information, the normal only two. For the octopole 
the complete information is contained in the three normals, 
but the dot-product statistic misses several degrees of free- 
dom.) We have estimated that including this extra struc- 
ture strengthens the 99% C.L. bound on the correlation of 
the quadrupole-octopole with the ecliptic plane to between 
99.93% C.L. (for the ILC map) and 99.996% C.L. (for the 
LILC map). 

In this work, we have continued to use cleaned, full- 
sky maps produced from the first year WMAP data. The 
concern with using these maps is the potential for residual 
Galactic foregrounds biasing the results. Though it is diffi- 
cult to see how Galactic contaminations can lead to eclip- 
tic correlations, we have studied this in two different ways: 
we have explored the properties of the foreground multipole 
vectors (in Secinj and we have explored how the multipole 
vectors and our results change when we perform a symmet- 
ric cut across the Galactic (and other) planes (in Sec. |7|l. 
As expected a Galactic foreground is dominated by the I20 
and Re(y3i ) modes. These are very different from the modes 
that dominate in the full-sky maps. As shown in Fig. |S|this 
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corresponds to the multipole vectors and normals for the 
full-sky maps being in very different locations than those 
for the foregrounds. As foreground contamination is slowly 
added to the full-sky maps we see how the full-sky multipole 
vectors move to the foreground multipole vectors. We have 
found that large foreground contaminations (jcj « 0.3 for 
the quadrupole, |c| « 1-3 for the octopole) are required to 
make the full-sky multipole vectors look like those from the 
known foreground maps. 

An alternative to using foreground maps is to mask out 
all information in the regions of the sky dominated by fore- 
grounds. This is a more conservative approach but throws 
away information about the CMB in some regions of the 
sky. For this reason any results from such a cut sky analy- 
sis will be weaker than the corresponding full-sky analysis. 
In Sec. we considered symmetric cuts across the Galac- 
tic (and other) planes to access the effect on the multipole 
vectors and correlations we have found from the full-sky 
analyses. As seen in Fig. |7] even small cuts lead to large 
uncertainties in our results. This is true independent of the 
plane about which we cut. The correlations we report from 
the full sky remain consistent in the cut sky analysis but 
are weakened as expected. This result is consis tent w i th the 
power equalization rec onstruction byfe iclcwicz. et al.1 ll2005h 
the cut sky analysis of ISlosar fc SeliakI 1120041) . 

As a final comparison of the quadrupole and octopole 
alignments we calculated the correlations for COBE maps 
(see Sec.lHJ. A gain the results are consistent with those from 
the WMAP full-sky maps (as shown in Fig. |HJ but do not 
have the statistical significance. 

The lack of power on the largest angular scales first ob- 
served by COBE and more recently confirmed by WMAP 
has motivated much of the study of the low multipole mo- 
ments, in particular the quadrupole and octopole. Ifowever, 
these are not the only multipole bands where there are pe- 
culiar features in the power spectrum. We have extended 
some of the studies to higher multipoles (Sees. 1^31 and IF!^ . 
Our tests suggest that there may exist peculiarities in these 
multipole ranges not solely in the power, but also in the 
structure of the multipoles. These studies, however, are not 
complete and thus it is not possible to assign statistical sig- 
nificance to them. They do, however, point the way for fu- 
ture work. 

To conclude, using the multipole vector decomposition 
we have shown that the quadrupole and octopole of the mi- 
crowave background sky are correlated with each other at a 
level that is excluded from being chance in excess of 99%. 
This comes about from a preponderance of peculiar correla- 
tions and is statistically independent of their observed lack 
of power. This observation is in bold contradiction to the 
predictions of pre-existing cosmological model, and argues 
against an infiationary origin for these fluctuations. In ad- 
dition, there is strong evidence (again of greater than 99% 
confidence) that the microwave background at these multi- 
poles is correlated with the geometry and direction of motion 
of the solar system. The observed signal is most unlikely to 
be due to residual contamination of the full-sky microwave 
background maps by known Galactic foregrounds. 

These results strongly suggest that either there is ad- 
ditional, unexplained foreground contamination of the mi- 
crowave background, potentially from a source local to our 
solar system or its neighborhood, or that there is an un- 



expected systematic error in the WMAP maps. We remain 
convinced by the WMAP team's ar guments that the re is 
no unexpected system atic error tBennett et aljl2003al : see 
also lFinkbeineill2004) . In particular, it is very hard to see 
how a north-south ecliptic asymmetry, or a quadrupole- 
octopole plane perpendicular to the ecliptic could be induced 
in the WMAP instrument or analysis pipeline. COBE, with 
largely independent error sources, saw compatible correla- 
tions. There is also the tantalizing suggestion by Archeops 
of a deficit in power near the first peak which is localized on 
the sky to the region of the ecliptic north pole. 

The astute reader will note that we have persisted in our 
failure to offer either a satisfactory possible explanation for 
an ecliptic-correlated foreground (especially one apparently 
concentrated in a plane perpendicular to the ecliptic) or a 
prediction that can be convincingly tested. Both are failings 
which we intend to remedy in the near future. However, we 
note that should indeed the low-£ microwave background 
prove to be dominated by a new foreground, this would im- 
ply that, barring an unexpected foreground alignment, the 
power in the underlying cosmic contribution at these mul- 
tipoles is likely to be suppressed below even the currently 
observed too-low value. It is at least amusing to note that 
the scale on which the lack of large scale correlations is then 
manifested is comparable to the horizon scale at the onset 
of cosmic acceleration. At the least this profound lack of 
large-angle correlations would further challenge generic in- 
fiationary models, maybe even general relativity on the scale 
of the observable universe will need to be reconsidered. 

Whatever the origin of these low-£ correlations, it 
is clearly necessary to reconsider any inferences drawn 
from the low-^ WMAP data, including the temperature- 
polarization cross-correlation. For example, our work sug- 
gests that the evidence for early reionization of the universe, 
resting as it does on the low-^ TE angular power spectrum 
should be viewed with a healthy dose of skepticism. 

On the experimental side, we are eagerly waiting for 
two new major and largely independent probes of the large- 
angle microwave background radiation that may shed new 
light on the anomalies discussed in this paper: polarization 
measurements by WMAP (though they are expected to be 
systematics-dominated on large scales), and meeisurements 
of temperature and polarization by the Planck experiment. 
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